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INTRODUCTION 

Of  recent  interest  in  breast  cancer  therapy  is  blocking  the  activity  of  HER2/neu  (a  member  of  the 
EGFR  family),  known  to  be  overexpressed  in  20-30%  of  breast  cancers  correlating  with  poor 
prognosis.  Another  major  concern  in  breast  cancer  is  the  status  of  p53-MDM2  where  MDM2  (a 
protein  known  to  downregulate  p53)  is  overexpressed  in  68%  of  malignant  breast  cancers  partly 
due  to  overexpression  of  the  estrogen  receptor.  Both  the  HER2/neu  and  MDM2  pathways  are 
regarded  separately,  but  it  is  possible  in  the  plethora  of  cellular  signaling  cascades  that  one  can 
have  an  effect  on  the  other  and  vice  versa.  It  is  evident  that  protein-protein  interactions  are 
important  in  all  signaling  cascades.  These  interactions  are  very  specific,  selective,  and 
evolutionarily  optimized.  Thus,  a  selective  interruption  of  protein-protein  interactions  by  small 
molecule  drugs  targeting  to  HER2/neu  and  MDM2  may  be  an  innovative  approach  to  cancer 
therapy.  We  are  developing  small-molecule  drugs  against  both  targets  and  hypothesize  that 
such  drugs  are  synergistic  when  used  in  combination. 
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BODY 

Our  main  task  was  to  employ  computational  methods  to  identify  small  molecules  able  to  fit  into 
the  hydrophobic  cleft  of  MDM2  and  to  evaluate  in  vitro  activity  of  these  compounds  in  p53  wild- 
type  and  mutated  cell  lines.  In  this  work,  we  performed  pharmacophore  model  studies  including 
the  substructure  concept  and  the  explicit  shape  constrains  to  the  active  conformation  of  Nutlin 
and  the  p53  fragment  that  binds  to  MDM2.  Database  screening  using  these  models  resulted  in  a 
series  of  structurally  novel  compounds.  Ninety  representative  hits  were  selected  for  initial 
cytotoxicity  assay  against  a  panel  of  breast  cancer  cell  lines  as  well  as  HCT 1 1 6  p53"^^"^  and  p53'^' 
cell  lines.  Eight  of  the  compounds  showed  desirable  activities  in  these  cell  lines. 

Functional  Feature  Pharmacophore  Model.  Functional  feature  pharmacophore  model  studies 
were  extensively  performed  on  both  Nutlin  and  the  p53  side  chain  of  the  conserved  triad  motif, 
F19,  W23,  L26  by  Catalyst  software  package.  The  conformation  of  Nutlin  and  p53  was  taken 
from  the  x-ray  structures  and  modeled  conformation.  For  Nutlin  pharmacophore  model 
development,  we  first  mapped  H-bond  donor,  H-bond  acceptor,  hydrophobic  and  aromatic  ring 
features  to  the  active  conformation.  We  then  assigned  geometrical  constraints  to  each  feature, 
i.e.,  coordinate  and  size  of  the  feature.  Finally,  all  the  selected  features  were  merged  into  one 
pharmacophore  model.  Following  a  similar  procedure,  a  pharmacophore  model  of  the  p53  triad 
motif,  F19,  W23,  L26,  was  also  derived. 

Shape-merged  single  compound  pharmacophore  query.  Catalyst  software  package  was 
used  to  generate  the  explicit  shape  constraint  of  Nutlin  or  p53  side  chain  ligand  (Figures  1  and 
2).  First,  the  active  conformation  of  the  template  molecule  (Nutlin  and  p53  fragment  of  amino 
acid  residues  from  19-26)  was  mapped  to  the  corresponding  feature  model,  respectively.  Then 
the  explicit  shape  constraint  was  generated.  Finally,  the  shape  was  merged  with  the  feature 
model  as  a  single  query  in  Catalyst,  that  we  referred  to  as  the  shape-merged  pharmacophore 
model.  The  similarity  tolerance  of  shape  is  an  important  parameter  to  adjust  the  percentage  of 
hits  yielded  from  database  screening.  The  values  of  0.5-1  was  used  as  default  value  in  the 
shape-merged  p53  side  chain  model,  but  for  the  shape-merged  Nutlin  feature  model,  we 
relaxed  the  shape  similarity  to  values  of  0.45-1 . 

Docking  studies.  The  docking  studies  were  carried  out  by  the  GOLD  (version  1.2  Genetic 
Optimization  for  Ligand  Docking)  software  package  running  on  our  multi-processor  linux  PC  and 
a  24-processor  Silicon  Graphics  Onyx  workstation  as  described.  We  used  the  crystal  structure 
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of  Nutlin  bound  MDM2  complex  (pdb  code  1RV1)  as  the  template  target,  and  the  ligand,  Nutlin, 
was  subsequently  removed  to  keep  the  binding  pocket  available.  Hydrogen  atoms  were  added 
to  both  protein  and  ligand  according  to  the  protonation  state  at  pH  7.0  during  the  docking 
studies. 

The  docking  area  was  defined  as  a  sphere  with  the  radius  of  15  A  at  the  central  coordinate 
position  (47.47,  12.31,  35.24)  of  the  Nutlin  molecule  as  observed  in  the  crystal  structure  of 
MDM2  binding  pocket.  The  ligand  was  placed  within  the  active  site  using  a  least-squares  fitting 
procedure.  All  docking  runs  were  performed  using  standard  default  settings  with  a  population 
size  of  100,  a  maximum  number  of  100,000  operations,  a  mutation  and  crossover  rate  of  95. 
The  scoring  function  was  contributed  basically  from  H-bond  interaction,  van  der  Waals 
interactions  within  the  complex,  and  the  ligand  internal  energy,  which  was  summarized  by  ligand 
steric  and  torsional  energies.  Explicit  electrostatic  interaction  was  ignored  but  was  modeled  into 
the  H-bond  interactions. 

Nutlin-based  pharmacophore  model.  Referring  to  the  interactions  between  Nutlin  and  the 
target  protein,  MDM2,  determined  in  the  X-ray  structure  (1RV1),  only  hydrophobic  features  (PI¬ 
PS)  were  considered  in  the  pharmacophore  model  (Figure  la).  The  features  PI,  P2,  and  P4 
were  mapped  by  the  benzene  ring  and  P3  by  aliphatic  methane  group,  which  was  observed  to 
interact  with  MDM2  by  its  Van  der  Waals  contacts.  Additionally,  P5  and  P6  were  recognized  by 
the  bromine  substitutes.  We  generated  explicit  shape  of  Nutlin  in  its  active  orientation,  as  shown 
in  Figure  1b.  Apparently,  both  feature  model  (Figure  la)  and  shape  constraint  itself  (Figure  1b) 
can  be  used  independently  as  searching  query  to  screen  small  molecule  database.  In  this  work 
however,  we  combined  the  shape  constraint  onto  the  feature  model  to  generate  a  shape- 
merged  pharmacophore  model,  labeled  as  SMI  (Figure  1b),  in  order  to  comprehensively 
describe  the  chemical  characteristics  of  Nutlin.  Consequently,  the  shaped-merged  model,  SMI, 
was  applied  to  screen  our  in-house  database  and  705  compounds  were  identified  to  fulfill  this 
query.  We  selected  structurally  diverse  compounds  for  further  cell-based  studies. 
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Figure  1.  Pharmacophore  models  generated  from  active  conformation  of  Nutlin.  (a)  Mapping 
Nutlin  onto  the  feature  pharmacophore  model,  (b)  Mapping  Nutlin  to  the  shape-only  query,  (c) 
Mapping  Nutlin  onto  shape-merged  pharmacophore  model. 


p53  side  chain-based  pharmacophore  modei 

Since  the  p53-MDM2  fragment  was  resolved  by  x-ray  technique,  several  peptide  libraries  were 
previously  designed  to  compete  with  the  p53  binding.  The  first  published  putative  p53-based 
pharmacophore  model  addressed  the  hydrophobic  features  recognized  by  the  side  chain  of 
F19,  W23  and  L26,  and  H-bond  donor  mapped  by  the  W23  indole  side  chain.  Its  application  to 
screen  the  NCI  database  resulted  in  a  series  of  non-peptide  sulfonamide  inhibitors.  In  this  work, 
however,  we  extended  the  pharmacophore  studies  by  taking  into  account  explicit  shape  of  the 
active  conformation  of  the  p53  fragment,  determined  by  the  x-ray  structure.  Two  shape-merged 
pharmacophore  models  were  systematically  generated  as  shown  in  Figure  2.  Figure  2a  shows 
the  model  with  the  consideration  of  only  three  hydrophobic  features,  but  enhanced  by  shape 
constraint,  referred  to  as  SM2.  Figure  2b  shows  the  shape-merged  model,  referred  to  as  SM3, 
that  keeps  additional  hydrogen-bond  feature.  Therefore,  SM2  and  SM3  were  treated  as 
independent  queries  to  virtually  screen  our  in-house  database. 
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Figure  2.  Pharmacophore  models  derived  from  the  side  chains  of  p53  fragment  interacting  with 
MDM2  p53  binding  domain,  (a)  p53  fragment  (amino  acide  residue  19-26)  mapping  onto  the  3- 
hydrophobic  feature  pharmacophore  model,  (b)  p53  fragment  (amino  acide  residue  19-26) 
mapping  onto  the  hetero-feature  pharmacophore  model. 

Database  search.  Substructure  queries,  and  the  shape-merged  pharmacophore  models 
derived  from  either  the  Nutlin  or  the  side  chain  of  p53  fragment  in  their  active  conformation, 
respectively,  were  applied  to  screen  our  in-house  database  by  Catalyst.  For  each  model  we 
identified  anywhere  from  700  to  900  hits.  Testing  selected  compounds  from  two  of  the  searches 
generated  novel  compounds  with  significant  activities  in  two  different  cytotoxicity  assays. 
Among  the  tested  compounds,  MP15,  MP17,  MP22,  and  MP31  showed  superior  activity  against 
cells  with  wild-type  p53  expression.  The  IC50  values  range  from  1-10  pM  using  MTT  assay  and 
0.5-5  pM  using  colony  formation  assay.  The  docked  structures  of  selected  compounds  are 
shown  in  Figure  3. 
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Figure  3.  Predicted  binding  mode  anaiysis  of  active  compounds.  MDM2  binding  surface 
is  dispiayed  in  yeiiow  (dot  or  soiid)  surface,  and  iigand  is  rendered  as  stick  or  baii-stick. 
(a)  and  (b)  show  the  MP15  binding  into  the  MDM2  pocket,  (a  is  for  top  view  and  b  for  side 
view),  (c)  and  (d)  show  the  MP17  binding  into  the  MDM2  pocket,  (e)  and  (f)  show  the  MP22 
binding  into  the  MDM2  pocket,  (g)  and  (h)  show  the  MP31  binding  into  the  MDM2  pocket. 
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Novel  inhibitors  of  HER2/neu 

For  HER2/neu  part  of  the  project  we  identified  several  novel  compounds  by  screening  a  subset 
of  our  in-house  collection  of  1 ,000  structurally  diverse  compounds.  Four  compounds  Lnz  1 0, 

18,  25,  and  28  that  showed  significant  inhibition  of  HER2/neu  kinase  activity  were  tested  against 
SKBR3  cells  using  MTT  (not  shown)  and  colony  formation  assay  (Figures  4  and  5). 


LNZ28  Treatment  Induces  a  Go/Gi-Phase  Arrest.  Cell  cycle  perturbations  induced  by  Lnz  28 
were  examined  in  SKBR3  cells.  The  analysis  of  DNA  profiles  by  flow  cytometry  indicated  that 
Lnz  28  induced  Go/Gi-phase  arrest.  As  shown  in  Figure  6,  nearly  81%  of  the  cells  were  still 
retained  in  Go/Grphase  after  48  hours  of  treatment  with  Lnz  28  (5  pM).  Similar  effects  were 
obtained  on  asynchronus  MDA-MB-435  cancer  cell  lines  (data  not  shown). 


noaeiea events: 

Ai  I  cycle  events:  1 7848 
Cycle  events  per  chennel:  1 82 
RCS:  3.580 
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Contro 
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58% 

67% 
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19% 
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4% 
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23% 
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Figure  6.  Cell  cycle  profile  of  cells  treated  with  5  pM  of  Lnz  28. 


In  collaboration  with  Ya-Qiu  Long  we  have  prepared  500  mg  of  Lnz  28  for  our  proposed  animal 
studies.  These  studies  will  be  performed  in  2008  through  different  funding  mechanism. 

In  parallel  we  have  performed  computational  search  of  our  5,000,000  compounds  and  identified 
additional  novel  HER2  antagonist.  The  full  disclosure  of  this  study  is  attached  in  the  appendix. 
This  paper  was  communicated  to  the  Journal  of  Medicinal  Chemistry  and  was  recently  accepted 
with  minor  revision.  Please  see  appendix. 
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KEY  RESEARCH  ACCOMPLISHMENTS 

Three  highly  robust  pharmacophore  models  were  developed  based  on  p53  peptide  fragment, 
Nutlin  structure,  and  MDM2  pocket.  These  pharmacophore  models  were  validated  against  a 
database  of  all  patented  MDM2  antagonists.  This  allowed  us  to  write  the  most  comprehensive 
review  of  inhibitors  of  p53-MDM2  interaction,  which  was  recently  published  in  Expert  Opinions  in 
Therapeutics  Patents. 

Several  small-molecule  inhibitor  of  p53-MDM2  interaction  were  identified.  One  example 
belongs  to  an  entirely  novel  class  of  anticancer  drugs.  We  filed  a  provisional  patent  on  these 
molecules. 

To  the  best  of  our  knowledge,  the  application  of  shape  merged  with  pharmacohporic  feature  and 
its  use  to  MDM2  as  a  target  has  never  been  previously  attempted.  This  could  potentially  be  a 
very  useful  tool  for  other  studies  targeting  protein-protein  interaction. 

We  also  screened  a  subset  of  our  in-house  collection  of  small-molecule  compounds  against 
HER2/neu.  A  series  of  novel  compounds  were  identified  and  subsequently  tested  in  cell-based 
assays.  We  also  wrote  a  very  comprehensive  review  of  inhibitors  of  HER2/neu,  which  was 
recently  submitted  to  Expert  Opinions  in  Therapeutics  Patents. 

At  least  four  lead  molecules  targeting  HER2-neu  activity  were  identified.  Structure-optimization 
in  collaboration  with  Dr.  Ya-Qiu  Long  was  initiated  and  we  have  prepared  500  mg  of  our  most 
potent  compound  LNZ-28  for  in  vivo  studies. 

We  also  performed  more  robust  computational  studies  and  identified  a  series  of  novel  HER2 
antagonists.  Account  of  this  study  is  summarized  in  our  recent  paper,  which  has  been  accepted 
by  the  Journal  of  Medicinal  Chemistry  with  minor  revision. 


Page- 1 1  - 


BC045544 


REPORTABLE  OUTCOMES 


Patent  Application 

Neamati,  N.  and  J.  Deng,  Structure-based  anticancer  drug  design  targeting  p53-MDM2 
interactions  (2005). 

Peer-reviewed  publications 

Deng,  J.;  Dayam,  R.;  Neamati,  N.,  Patented  small  molecule  inhibitors  of  p53-MDM2  interaction. 
Expert  Opin  Ther  Patents  2006,  16,  (3),  165-188. 

Dayam  R,  Grande,  F.;  Al-Mawsawi,  L-Q.,  Neamati,  N.  Recent  advances  in  the  design  and 
discovery  of  small-molecule  therapeutics  targeting  HER2/neu.  Expert  Opin  Ther  Patents  2007, 
17,  (1),  83-102. 


Gundia,  R.;  kazemi,  R.;  Sanam,  R.;  Muttineni,  R.;  Dayam,  R.;  Jagarlapudi,  ARP,;  Neamati,  N. 
Discovery  of  novel  small-molecule  inhibitors  of  HER2/neu:  Combined  ligand-based  and  target- 
based  approach.  J.  Med.  Chem.  (Accepted  with  minor  revision). 


Deng,  J.;  Taheri,  L.;  and  Neamati,  N.  Design  and  discovery  of  novel  small-molecule  compounds 
targeting  p53-MDM2  interaction  (Under  preparation) 


Page- 12  - 


BC045544 


CONCLUSIONS 

We  have  successfully  discovered  a  series  of  novel  inhibitors  of  both  p53-MDM2  interaction  as 
well  as  HER2/neu  kinase  activity.  This  success  attests  to  your  unique  approach  in  structure- 
based  drug  design.  Our  compounds  show  significant  cytotoxicity  against  breast  cancer  cell 
lines.  Our  animal  studies  will  be  performed  in  the  2008  under  different  funding  mechanism. 
Thus  far,  we  have  filed  one  patent  application  and  published  three  papers  and  an  additional 
manuscript  is  under  preparation.  This  is  a  major  project  in  our  laboratory  and  our  goal  is  to  use 
these  preliminary  results  to  apply  for  a  bigger  grant. 
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Patented  small  molecule  inhibitors 
of  p53-MDM2  interaction 

Jinxia  Deng,  Raveendra  Dayam  &  Nouri  Neamati^ 
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1985  Zonal  Avenue,  PSC304A,  Los  Angeles,  California  90089,  USA 

The  interaction  between  p53  and  murine  double  minute  2  (MDM2)  provides^ 
an  attractive  drug  target  in  oncology.  Small  molecule  inhibitors  of  this  inter¬ 
action  have  not  only  provided  strong  evidence  for  blocking  the  protein-pro¬ 
tein  interaction,  but  are  also  extremely  useful  as  biological  probes  and 
ultimately  as  novel  therapeutics.  Here,  a  comprehensive  review  of  the  pat¬ 
ented  small  molecule  inhibitors  of  the  p53-MDM2  interaction  are  provided. 
These  inhibitors  are  divided  into  eight  classes  of  compounds  that  include 
c/s-imidazolines,  benzodiazepines,  fused  indoles,  substituted  piperazines, 
substituted  piperidines,  aryl  boronic  acids,  spiro-indoles,  and  a-helix  mimetic 
compounds.  The  best  documented  class  of  compounds,  c/s-imidazolines  (e.g., 
Nutlins)  are  selective  and  potent  inhibitors  of  the  p53-l\/IDM2  interaction, 
and  selected  examples  exhibit  potency  in  the  nanomolar  range.  Nutlins 
induce  apoptosis  in  p53  wild-type  cells  and  show  in  vivo  efficacy  in  mice 
xenograft  models.  Additional  strategies  briefly  discussed  in  this  review,  and 
which  are  under  current  exploration  in  targeting  the  p53  pathway,  include 
the  inhibition  of  MDM2-mediated  p53  ubiquitylation  and  restoration  of 
DNA-binding  activity  of  mutant  p53  protein  using  small  molecules. 
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.  Introduction 

- 

1.1  MDM2 

Murine  double  minute  2  (MDM2)  is  an  oncogene  originally  isolated  from  a  mouse 
fibroblast  cell  line  that  underwent  a  spontaneous  transformation  [i].  Later  this  gene 
was  rediscovered  as  a  p53  binding  protein  in  several  rat  fibroblast  cell  lines  [2].  The 
human  homologue  has  since  shown  to  be  overexpressed  in  8%  of  various  sarcomas, 
up  to  68%  in  angiosarcomas  and  ~  30%  in  soft  tissue  sarcomas  [3-7].  A  large  number 
of  cancers  with  MDM2  overexpression  harbour  wild-type  p53.  In  some  cancers, 
MDM2  overexpression  can  be  correlated  with  accelerated  cancer  progression  and 
poor  response  rate  [8-i3].  Several  studies  suggest  that  MDM2  overexpression  is  associ¬ 
ated  with  increased  metastasis,  although  the  function  and  mechanism  of  MDM2  in 
cancer  metastasis  demands  more  studies  [io,i4,i5].  Three-dimensional  structural  infor¬ 
mation  is  not  yet  available  for  the  full-length  MDM2  protein.  The  domain  composi¬ 
tions  of  p53  and  MDM2  are  shown  in  Figure  1.  However,  biological  mapping  has 
shown  that  the  human  MDM2  gene  encodes  for  a  protein  containing  491  amino 
acids  with  multiple  domains  and  consensus  sequences.  These  include  a  p53-binding 
domain,  a  nuclear  localisation  signal  (NLS),  a  nuclear  export  signal  (NES),  an  acidic 
region,  a  zinc-finger  and  a  ring-finger  domain  (Figure  lA)  [6,16,17].  MDM2  is  largely 
considered  a  major  negative  regulator  of  p53  tumour  suppressor  gene  [I8,i9]. 

1.2  p53 

p53  contains  393  amino  acid  residues.  The  protein  is  a  transcription  factor  that  regu¬ 
lates  the  expression  of  proteins  that  mediate  distinct  cellular  responses  such  as 
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Figure  1.  Schematic  representations  of  the  functional  domains  of  p53  and  MDM2.  A)  Domain  composition  of  p53  and  MDM2. 
B)  Core  domain  of  p53  complexed  with  the  DNA  (PDB  1TSR  [26]).  Image  generated  by  Molsoft  ICM-Browser.  C)  Tetramerisation  domain 
of  p53  (PDB  1C26  [31]).  Image  generated  by  Molsoft  ICM-Browser.  D)  p53  peptide  bound  MDM2  (PDB  1YCR  [35]).  Green  surface 
represents  MDM2.  Helix  in  red  represents  p53  peptide  bound  to  MDM2,  and  the  side  chains  in  orange  show  triad  motif:  L26,  W23,  FI  9. 
Image  generated  by  PyMOL.  Colour  figure  available  online  at  www.ashley-pub.com. 

NES:  Nuclear  export  signal;  NLS:  Nuclear  localisation  signal. 


apoptosis  and  cell  growth  arrest  120,21].  Not  surprisingly,  p53  is 
the  most  frequently  inactivated  protein  in  human  cancers  [22], 
Over  50%  of  human  cancers  contain  p53  mutations  that  inac¬ 
tivate  its  transcription  regulation  activity  [23].  Most  cancer- 
associated  p53  mutations  occur  within  the  DNA-binding 
domain,  rendering  p53  unable  to  act  as  a  sequence-specific 
DNA  jawing  transcription  factor  [24].  The  fact  that  muta- 
tionavractivation  of  p53  is  observed  as  one  of  the  most  com- 
cancer  genetic  events  indicates  that  p53  plays  a  critical 
Q^e  in  suppressing  the  tumour  growth.  Structural  studies  have 
’  shown  that  p53  is  made  up  of  at  least  four  functional  domains 
that  regulate  its  transcriptional  activities:  i)  an  N-terminal 
transactivation  domain,  which  is  required  for  the  interaction 
with  transcriptional  protein  machinery  as  well  as  a  proline-rich 
domain;  ii)  a  central  conserved  DNA-binding  core  domain;  iii) 
a  tetramerisation  domain  assisting  in  sequence-specific  DNA 
binding;  and  iv)  a  C-terminal  negative  regulatory  domain, 
which  when  phosphorylated,  primes  the  latent  sequence-spe¬ 
cific  DNA-binding  function  of  p53  for  activation  [6].  The 
N-terminal  transactivation  domain  contains  the  binding 


pocket  where  most  of  the  important  negative  cellular  regula¬ 
tors  of  p53  bind,  that  is  MDM2  oncoprotein  (Figure  lA)  [25]. 
The  structures  of  the  DNA  binding,  tetramerisation  and  the 
C-terminal  regulatory  domains  of  human  p53  have  been  deter¬ 
mined  by  NMR  or  crystallography  (Figures  IB  and  C)  [26-34]. 
Currently,  there  is  no  such  structural  information  available  on 
the  entire  human  p53  N-terminal  domain.  However,  the  crys¬ 
tal  structure  of  a  p53  peptide  that  contains  residues  15  —  29  in 
complex  with  MDM2  has  been  resolved  (Figure  ID)  [35]. 

1.3  P53-MDM2  interaction 

MDM2  regulates  p53  protein  activity  in  at  least  three  differ¬ 
ent  ways  [36].  Briefly,  i)  it  binds  to  the  p53  transactivation 
domain,  thereby  blocking  the  transcriptional  activity;  ii)  sub¬ 
sequent  to  p53-binding,  it  induces  nuclear  export;  and  iii)  it 
stimulates  the  degradation  of  p53  by  catalysing  ubiquitylation 
of  the  protein  through  the  ubiquitin/proteasome  pathway  (for 
recent  reviews,  see  [37,38]).  Transcription  of  the  MDM2  gene 
itself  is  regulated  by  p53  in  response  to  cellular  stress,  indicat¬ 
ing  that  p53  and  MDM2  form  an  autoregulatory  feedback 
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Figure  2.  Simplified  model  of  autoregulation  between  p53  and  MDM2.  p53  mediates  the  l\/IDM2  expression,  and  in  turn  MDI\/I2 
regulates  p53  activity  through  blocking  its  transcriptional  activity,  promoting  its  nuclear  export,  and  stimulating  its  degradation  in  both 
nucleus  and  cytoplasm.  DNA  model  was  generated  by  ICM-Browser  software. 


loop  (Figure  2)  [36] .  Therefore,  as  one  of  the  approaches  to  sta¬ 
bilise  p53  protein  levels,  MDM2  has  been  selected  as  a  very 
attractive  target  for  cancer  therapeutics.  Efforts  are  underway 
to  design  peptides  or  small  molecule  compounds  that  disrupt 
the  p53-MDM2  interaction.  Biological  and  functional  stud¬ 
ies  have  mapped  residues  19  —  102  ofMDM2  as  the  domain 
for  p53  binding,  while  the  MDM2  binding  domain  of  p53 
was  narrowed  down  to  a  15 -amino  acid  residue  peptide 
(amino  acid  residues  15  -  29;  Figure  ID)  [6].  The  structure 
determination  of  MDM2  with  the  15-mer  wild-type  p53  pro¬ 
vided  a  precise  description  of  the  p53-MDM2  interaction. 
This  revealed  that  MDM2  possesses  a  well-defined  deep 
hydrophobic  binding  pocket.  A  triad  of  amino  acid  residues 
FI 9,  W23  and  L26  from  MDM2  binding  15-mer  p53  pep¬ 
tide  occupies  the  deep  hydrophobic  pocket  (Figure  ID)  [35].  A 
key  H-bond  interaction  between  the  indole  nitrogen  of  W23 
in  p53  and  the  target  MDM2,  and  several  van  der  Walls  con¬ 
tacts  involved  in  the  stabilisation  of  the  complex  were 
observed.  The  structure  of  p53-bound  MDM2  provided  a 
rational  template  for  structure-based  design  of  small  molecule 
inhibitors.  In  the  current  review,  the  first  focus  will  be  on  the 
patented  small  molecule  inhibitors  of  the  p53— MDM2 


interaction.  Second,  recent  progress  towards  the  discovery  of 
small  molecule  inhibitors  that  prevent  MDM2-mediated  p53 
ubiquitylation  will  be  covered  briefly.  Finally,  the  discovery  of 
small  molecules  that  have  the  ability  to  reactivate  the  DNA 
binding  and  apoptosis-inducing  function  of  mutant  p53 
protein  will  be  described. 

2.  Discovery  of  small  molecule  inhibitors  of 
p53-MDM2  interaction 

Most  of  the  early  studies  employed  peptides  rather  than  small 
molecules  to  disrupt  the  p53-MDM2  interaction.  These  initial 
studies  suggested  that  the  peptides  derived  from  p53  could  be 
used  as  probes  to  investigate  this  specific  interaction.  They  also 
confirmed  the  critical  function  of  the  triad  motif  as  revealed 
from  the  p53-MDM2  complex  [39-4i].  Subsequently,  a  potent 
12-mer  peptide  was  identified  by  screening  phage  peptide 
libraries  [42,43].  Other  shorter  peptides  showed  tighter  binding 
as  observed  from  the  kinetic  and  thermodynamic  studies  [44]. 

Several  small  molecule  inhibitors  of  the  p53— MDM2  inter¬ 
action  that  include  chalcones  and  chlorofusin  have  previously 
been  discovered  [45-50].  However,  major  developments  in  the 
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field  were  the  discoveries  of  Nutlins  (representatives  of  cw-imi- 
dazolines)  and  benzodiazepine  compounds  [50,51].  The  subse¬ 
quent  co-crystal  structures  of  these  compounds  with  MDM2 
greatly  helped  in  the  understanding  of  critical  binding  inter¬ 
actions  between  small  molecules  and  the  MDM2  protein. 
These  studies  established  that  protein-protein  interactions 
can  be  successfully  manipulated  by  small  molecule  inhibitors. 
Therefore,  the  p53  pathway  can  be  effectively  studied  at  the 
molecular  level  to  develop  novel  anticancer  drugs  selectively 
targeting  its  critical  interaction  with  MDM2. 

Several  reviews  have  summarised  the  early  progress  in  the 
design  of  peptides  and  small  molecule  inhibitors  targeting 
p53-MDM2  interaction  [19,43,52-55].  However,  the  emergence 
of  patented  small  molecules  in  the  past  5  years  has  not  been 
covered.  Herein,  the  entire  patent  literature,  published  since 
2000,  for  the  p53— MDM2  interaction  inhibitors  has  compre¬ 
hensively  reviewed.  It  is  notable  to  bear  in  mind  that  many 
compounds  claimed  to  be  inhibitors  of  p53— MDM2  interac¬ 
tion  have  yet  to  be  independently  validated.  These  patents 
were  not  discarded.  Additionally,  in  several  patents,  many 
compounds  were  disclosed  but  their  biological  activities  were 
not  presented.  Even  though  it  is  very  difficult  to  present  a 
coherent  structure-activity  relationship  (SAR)  at  this  time,  it 
was  decided  that  it  was  still  informative  to  present  selected 
representatives  from  each  study.  The  purpose  is  to  provide  a 
rapid  access  to  compound  structures  that  are  claimed  to  be 
inhibitors  of  p53— MDM2  interaction,  beneficial  to  interested 
researchers.  > 

In  general,  the  small  molecule  inhibitors  fall  into  eight 
structural  categories.  They  are  aV-imidazolines  (as  represented 
by  Nutlins  [50[),  benzodiazepines,  fused  indoles,  substituted 
piperazines,  substituted  piperidines,  aryl  boronic  acids,  spiro- 
oxindoles  and  a-helix  mimetic  compounds.  To  the  authors’ 
knowledge,  Nutlins,  examples  of  «f-imidazolines,  are  the  only 
selective  and  potent  inhibitors  that  have  been  studied  in  detail 
in  cell  culture  and  mice  xenograft.  Therefore,  the  review 
begins  with  the  discussior^^w-imidazolines  and  the  recent 
progress  in  the  design  .wMiJmscovery  of  potential  inhibitors  of 
p53— MDM2  intei^^^. 

2.1  C/s-imidaz(^tnes 

The  first  biological  application  of  a'i-imidazolines  could  be 
illustrated  by  TA-383  (cw-2-(4-chlorophenyl)-4,5-diphenyl-2- 
imjjazoline;  1,  Figure  3),  which  was  recognised  as  a  novel  anti- 
Cy^umatic  agent.  Later  the  compound  was  shown  to  enhance 
macrophage  migration  [56,57].  In  addition,  4,5-diarylimidazo- 
lines  were  recently  evaluated  as  inhibitors  of  the  P2Xy  recep¬ 
tor,  which  is  a  ligand-gated  ion  channel  that  is  present  on  a 
variety  of  cell  types,  including  mast  cells,  macrophages  and 
lymphocytes,  all  thought  to  be  involved  in  inflammation  and 
autoimmune  processes  [58]. 

Nutlins  are  the  representatives  of  a  class  of  af-2,4, 5-tri- 
phenyl-imidazolines,  and  it  is  only  in  recent  years  that  they 
were  recognised  as  the  first  potent  and  selective  small  mole¬ 
cule  inhibitors  of  the  p53— MDM2  interaction  [50,59,60,i0i-i04]. 


Cl 


Figure  3.  Structure  of  TA-383,  a  c/s-imidazoline  analogue-O 

These  compounds  were  identified  from  library  ss^^^^ng  of  a 
large  diverse  synthetic  chemicals  database  [50[ .  Three  reported 
Nutlins  showed  potency  for  blocking  the  p53— MDM2  com¬ 
plex  in  the  range  of  100  to  300  nM  with  ~  150-  to  200-fold 
difference  between  the  enantiomers  [50]. 

Nutlins  were  initially  examined  at  the  cellular  level  for  their 
effect  on  p53,  MDM2,  and  p21  protein  expression.  Results 
confirmed  that  Nutlin  treatment  accumulated  wild-type  p53 
levels  and  consistently  lead  to  an  elevation  ofMDM2  and  p21 
proteins  [50].  In  addition,  it  was  shown  that  the  increased  p53 
level  was  due  to  a  post-translational  mechanism  rather  than 
the  elevated  expression  of  the  p53  gene.  Cell  cycle  analysis  of 
bromodeoxyuridine  (BrdU) -labelled  cancer  cells  demon¬ 
strated  that  Gj  and  G2/M  phase  fractions  were  increased, 
whereas  cells  progressing  into  S-phase  were  completely 
deleted  after  one-day  Nutlin  treatment.  Results  from  a  3-(4,5- 
dimethylthiazol-2-yl)  -2,5  -diphenyltetrazolium  bromide 

(MTT)  assay  indicated  that  cytotoxic  activity  of  the  com¬ 
pounds  on  the  wild-type  p53  cells  was  significantly  higher 
than  mutant  p53  cells,  suggesting  that  the  p53  pathway  is 
only  turned  on  in  wild-type  p53-containing  cells.  The  active 
enantiomer  of  Nutlin  induced  apoptosis  in  50%  of  cells, 
whereas  the  inactive  enantiomer  had  no  effect.  Finally,  one 
analogue  was  investigated  in  osteosarcoma  SJSA-1  tumour 
xenografts  in  nude  mice.  The  20-day  oral  administration  of 
this  compound  effectively  inhibited  90%  tumour  growth  as 
compared  to  the  vehicle  control  [50,59]. 

To  investigate  the  Nutlin  binding  mode,  the  crystal  struc¬ 
ture  of  the  human  MDM2— Nutlin  complex  was  determined 
at  2.3  A  [50] .  Figure  4  shows  the  co-crystal  structure  of  a  Nutlin 
analogue  bound  to  MDM2  (PDB  IRVl).  It  clearly  shows 
that  the  Nutlin  analogue  projects  its  functional  groups  into 
the  p53  binding  pocket  of  MDM2  by  mimicking  the  confor¬ 
mation  of  the  p53  triad  motif,  that  is,  the  side  chain  of  FI 9, 
W23  and  L26. 

Some  representative  molecules  disclosed  from  recent  pat¬ 
ents  are  shown  in  Figures  5  and  6  [ioi-io4[.  N-substituted  imi¬ 
dazolines  are  shown  in  Figure  5,  whilst  un-substituted 
imidazolines  are  shown  in  Figure  6.  In  cell-free  and  cell-based 
assays,  these  molecules  were  observed  to  inhibit  the  p53- 
MDM2  interaction  with  a  potency  of  ~  100-fold  greater  than 
that  of  a  p53-derived  peptide.  In  addition,  the  cell-based 
assays  demonstrated  mechanistic  activity.  Incubation  of 
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Figure  4.  Crystal  structure  of  the  Nutlin  analogue  bound  MDM2  (PDB  1RV1  [50])i^lue  surface  represents  the  p53  binding  pocket  of 
MDI\/I2  (image  generated  by  PyMOL).  The  stick  model  shows  the  Nutlin  analogue.  Colour  figure  available  online  at  www.ashley-pub.com. 

cancer  cells  containing  wild-type  p53  resulted  in  p53  protein  .  receptor  for  HIV,  thereby  inhibiting  HIV  infection  [109].  Pre¬ 
accumulation,  p53-regulated  p21  gene  induction  and  cell  ^Viously,  a  class  of  4, 5, 6,7-tetrahydro-5-methylimidazo[4, 5,1- 


cycle  arrest  in  Gl,  G2  phase,  leading  to  potent  antiprolifer^ 
tive  activity  against  wild-type  p53  cells  in  vitro.  In  contrast, 
similar  results  were  not  observed  in  p53  mutant  cancer  cells  at 
comparable  drug  concentrations.  _  ^ 

Another  set  of  compounds  bearing  the^h-imdazoline 
group  were  discovered  as  inhibitors  of  the  p53-MDM2  inter¬ 
action  and  showed  antiproliferative  activities  [los-ios].  Several 
compounds  exemplified  by  I  and  II  (Tables  1  and  2)  were 
claimed,  and  the  representative  molecules  and  their  inhibitory 
activities  are  given  in  Tables  1  and  2.  From  a  simpler  formula, 
a  group  of  derivatives  containing  4, 5-dihydro- 1,4-dimethyl- 
2, 5-diphenyl- l//-imidazole  were  also  observed  to  induce  anti¬ 
proliferative  activity  through  inhibition  of  the  p53— MDM2 
interaction.  The  inhibitory  activity  of  each  compound  was 
measured  as  a  percentage  of  the  bound  MDM2  in  treated  ver¬ 
sus  untreated  wells.  The  IG50  values  of  the  disclosed  molecules 
A  from  the  invention  ranged  0.5  -  150  pM.  In  summary,  the 
discovery  of  Nutlins  as  potent  inhibitors  for  the  disruption  of 
p53  with  MDM2  in  vitro  and  in  vivo  provided  compelling 
evidence  that  such  molecules  have  excellent  therapeutic 
potential.  Therefore,  validation  of  p53-MDM2  disruption  as 
a  viable  therapeutic  target  provides  strong  support  for  target¬ 
ing  protein— protein  interactions  in  general.  This  is  a  clear  step 
forward  towards  novel  therapeutics  in  the  future. 

2.2  Benzodiazepines 

The  biological  activities  of  benzodiazepine  derivatives  are  very 
well  documented.  For  example,  some  analogues  have  been 
recognised  to  bind  the  chemokine  receptor  CXCR4,  a 


7^]benzodiazepin-2-(l//)-one  compounds  were  reported  pos¬ 
sessing  significant  inhibitory  activity  against  HIV-1  with  IGjq 
values  ranging  0.0034  -  238  pM  [61].  Separate  studies  showed 
tricyclic  diazepines  have  antiproliferative  activity  on  leukae¬ 
mia  cells  [110].  Donepezil  analogues  have  also  been  shown  as 
neurocyte  protective  agents  [111].  Recently,  the  imine-amine 
benzodiazepine  dimers  were  also  observed  to  efficiently  bind 
to  DNA,  and  exhibited  in  vitro  antitumour  activity  in  many 
cancer  cell  lines  that  included  leukaemia,  lung  cancer,  renal 
cancer  and  colon  cancer  [62-64].  Another  study  demonstrated 
that  the  benzothiadiazine-4-one  derivatives  were  effective 
anticoagulants  through  inhibition  of  Factor  Xa  activity  [65,66]. 
In  addition,  several  benzodiazepines  also  showed  anxiolytic 
and  herbicidal  activities  [67,68]. 

Although  the  biological  studies  of  benzodiazepine  deriva¬ 
tives  can  be  traced  to  the  early  1960s,  their  potential  as  anti¬ 
cancer  agents  has  only  recently  been  discovered.  Several  of 
these  compounds  were  found  to  effectively  block  the  p53- 
MDM2  interaction.  These  compounds  were  identified  by 
screening  a  synthetic  library  containing  >  22,000  compounds 
using  a  high-throughput  binding  assay  sensitive  for  MDM2 
binding  affinity  [69]. 

A  total  of  217  compounds  were  disclosed  in  two  recent  pat¬ 
ents  each  containing  a  critical  motif;  some  selected  com¬ 
pounds  are  shown  in  Figure  7  [112,113].  A  fluorescein  peptide- 
based  assay  was  used  to  measure  drug  potency.  Structures  of 
17  representative  compounds  are  shown  in  Figure  7.  Eleven 
compounds  inhibited  p53-MDM2  interaction  with  ICjq 
values  of  <  10  pM  (Table  3). 
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Figure  5.  Representative  examples  of  c/s-imidazoline  inhibitors  of  p53-l\/IDI\/l2  interaction. 
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Figure  6.  Representative  c/s-imidazoline  inhibitors  of  p53-l\/IDI\/l2  interaction. 
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A  modified  1 ,4-diazepine  scaffold  and  some  selected  deriv¬ 
atives  are  shown  in  Figure  8  [Ii4,ii5].  A  total  of  1 14  examples 
was  synthesised  and  tested  [ii3].  The  potency  of  these  com¬ 
pounds  ranged  from  0.05  to  >  100  pM.  The  activities  of  the 
selected  representative  compounds  are  given  in  Table  4. 

To  further  understand  the  molecular  detail  of  such  types 
of  benzodiazepinediones  binding  to  MDM2,  two  represent¬ 
ative  compounds  were  selected  for  X-ray  studies  [ii6].  The 
chemical  structures  of  the  two  compounds  are  shown  in 
Figure  9.  Coordinates  of  benzodiazepine  analogue  64  bound 
to  MDM2  have  been  determined  at  2.7  A  by  X-ray 


crystallography  (PDB  1T4E),  as  shown  in  Figure  10  [5i]. 
From  the  binding  conformations,  it  is  clear  that  both  Nutlin 
(see  Figure  3)  and  the  selected  benzodiazepine  analogue  64 
stay  in  a  similar  orientation  and  interact  with  MDM2  in  an 
orientation  that  mimics  the  side  chain  of  p53.  As  a  conse¬ 
quence,  the  p53  is  released  and  stabilised  in  the  cell.  Cell 
growth  suppression  was  measured  in  JAR  choriocarcinoma 
cells  that  overexpress  p53  and  MDM2.  Results  with  benzo¬ 
diazepine  analogue  64  treatment  confirmed  the  antiprolifer¬ 
ative  activities  were  directly  due  to  the  selective  interruption 
of  the  p53-MDM2  interaction.  In  addition,  the  discovery 
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Figure  7.  Representative  examples  of  benzodiazepine  inhibitors  of  p53-MDM2  interaction. 
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Figure  7.  Representative  examples  of  benzodiazepine  inhibitors  of  p53-MDM2  interaction  [Continued). 


Table  3.  Inhibition  of  MDM2  binding  to  p53  by 
representative  benzodiazepines. 


o'J' 


Compound 


Inhibition  of  MDM2  binding 
to  p53  (ICsoPM) 


42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 


1.7 

1.8 
13 
56 

7.8 
10.1 
11.7 
10.5 

1.9 

10 


O 

0.87 

>^1.70 


0.62 

0.47 

0.62 

0.58 


of  benzodiazepines  as  inhibitors  of  p53— MDM2  interaction 
was  the  first  report  of  such  types  of  compound  acting  as 
a-helix  mimetics  [5i].  Therefore,  this  finding  further  sug¬ 
gests  that  the  exploration  of  the  benzodiazepinedione  scaf¬ 
fold  towards  the  design  of  inhibitors  for  the  disruption  of 
protein-protein  interaction  is  a  valid  approach  [5i]. 

2.3  Substituted  piperazines 

Numerous  studies  on  the  biological  activities  of  4-phenyl- 
piperazine- 1 -carboxamide  substituents  have  been  reported. 
Some  have  been  shown  to  be  orally-active  nonsteroidal  andro¬ 
gen  receptor  inhibitors,  whilst  some  have  shown  antibacterial 
activity  [70-72].  A  new  synthetic  piperazine  derivative,  which 


inhibits  microtubule  activity,  k  reported  to  exhibit  anti- 
angiogenic  and  anticancer  acftVfies  [731.  Another  study  indi¬ 
cated  that  several  piperazine-based  synthetic  molecules  were 
potent  inhibitors  of  plasminogen  activator  inhibitor- 1,  a 
member  of  the  serine  protease  inhibitor  (serpin)  family  acting 
as  the  major  negative  regulator  of  tPA  and  uPA  [74],  Some 
piperazine^derivatives  were  recently  reported  as  histone 
deacetylase  and  NKl  inhibitors  [117-119]. 

Recently,  the  piperazine  derivatives  were  proposed  as  use- 
ful^^probes  of  the  p53-MDM2  interaction  and  as  novel 
z^agents  for  the  treatment  of  cancer  [120].  The  activities  were 
^  measured  by  a  modified  ELISA  using  a  histidine-tagged 
MDM2,  a  GST-p53  fusion  protein  and  a  nickel  chelate- 
alkaline  phosphatase  (NTA-AP). 

4-Phenylpiperazine  derivatives  are  another  example  of 
p53— MDM2  interaction  inhibitors  that  have  been  reviewed 
recently  [75,76,121].  The  IC50  value  of  these  compounds  ranged 
0.03  -  200  pM.  Compound  70  has  an  IC50  value  of  ~  4  pM. 
Representative  compounds  are  shown  in  Figure  1 1 .  These 
compounds  are  claimed  to  be  useful  in  combination  with 
other  anticancer  drugs  to  treat  tumours  where  MDM2  is 
overexpressed  [75]. 

2.4  Substituted  piperidines 

A  separate  discovery  provided  a  novel  series  of  piperidines 
with  cw-3,4-dialkoxy  substitutions.  These  compounds  were 
claimed  as  small  molecule  inhibitors  of  the  p53— MDM2 
interaction.  The  3,4-bis(benzyloxy)  piperidines  displayed  a 
variety  of  biological  activities.  For  example,  several  deriva¬ 
tives  were  useful  to  treat  Gaucher  disease,  a  genetic  disorder, 
and  cancer,  whereas  some  other  examples  were  claimed  to  be 
effective  in  the  prevention  and  treatment  of  Alzheimer’s  dis¬ 
ease  [122-124].  Moreover,  the  piperidine  motif  has  also  been 
recognised  as  a  core  structure  in  several  HIV  and  feline 
immunodeficiency  virus  (FIV)  protease  inhibitors  [125].  Some 
examples  were  reported  to  target  human  coronavirus  severe 
acute  respiratory  syndrome. 

Amongst  the  1 5  patents  covering  piperidines,  two  describe 
their  potential  as  inhibitors  of  the  p53— MDM2  interaction 
[126-127].  Representative  examples  are  shown  in  Figure  12.  The 
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Figure  8.  Representative  diazepine  inhibitflhrbf  p53-l\/IDI\/l2  interaction 

Table  4.  Inhibition  of  MDM2,bnY^g  to  p53  by 
representative  compound^  > 

- - T - 

Inhibition  of  l\/IDM2  binding 

to  p53  (IC50  pM) 


Compound 


”  dS’ 

- 

0.1  -  1.0 

q 

rn 

1 

q 

0.1  -  1.0 

0.1  -  1.0 

S 

q 

1 

q 

binding  orientation  of  the  selected  compounds  to  the  p53- 
binding  pocket  of  MDM2  was  validated  by  NMR  studies. 

2.5  Aryl  boronic  acids 

Several  aryl  and  heteroaryl  boronic  acids  were  recently  tested 
in  vitro  for  their  hormone-sensitive  lipase  inhibitory  properties. 
The  most  potent  compound  had  an  ICjq  value  of  17  nM  [77]. 
Another  study  showed  that  some  boronic  acids  were  potent 


inhibitors  of  EGFR  and  VEGFR-1  tyrosine  kinases  [78].  Other 
activities,  such  as  the  inhibition  of  P-lactamases  and  serine 
protease,  were  also  documented  [128-130]. 

A  recent  study  described  their  anticancer  activities  [131]. 
The  compounds  were  claimed  to  inhibit  MDM2  expression 
and  disrupt  the  p53-MDM2  protein  complex.  Structures  of 
four  selected  aryl  boronic  acids  are  shown  in  Figure  13.  The 
anticancer  activities  of  a  representative  compound  88  are 
given  in  Table  5. 

2.6  Fused  indoles 

A  series  of  fused  indoles  shown  in  Figure  14  were  reported  to 
interrupt  the  wild-type  p53— MDM2  interaction  [132].  These 
were  demonstrated  to  induce  cell  death  in  different  types  of 
human  tumour  cells,  for  example,  lung,  colon  and  breast  can¬ 
cer.  Moreover,  the  compounds  were  shown  to  stabilise  the  p53 
expression  and  prevent  its  degradation  in  living  cells.  Restora¬ 
tion  of  p53  transcriptional  activity  was  also  validated  by  the 
standard  biological  assay.  Selective  growth  suppression  activi¬ 
ties  of  a  set  of  representative  compounds  in  wild-type  p53- 
expressing  colon  cancer  cells  are  given  in  Table  6. 


176 


Expert  Opin.  Then  Patents  (2006)  16(2) 


Deng,  Dayam  &  Neamati 


64 


64a 


Figure  9.  Compounds  64  and  64a  were  co-crystallised  with  MDM2. 


qJttJFe  of  the  MDI\/I2  complexed  with 
[^ue  64  {PDB1T4E  [51]).  Blue  surface 


Figure  10.  Crystal  sti 
benzodiazepine 

represents  the  d^^Binding  pocket  of  l\/lDM2.  The  stick  model 
shows  the  b^m«iazepine  analogue  (image  generated  by  PyMOL). 
Colour  ^pre available  online  at  www.ashley-pub.com. 

j  2.7  Spiro-oxindoles 

^  From  the  biological  mapping  studies  and  X-ray  structural 
validation,  it  was  envisaged  that  W23  of  p53  appears  to  be 
most  critical  for  MDM2  by  forming  hydrophobic  and 
Id-bond  interactions.  A  recent  study  clearly  showed  that 
oxindoles  can  mimic  the  side  chain  of  W23  for  interaction 
with  MDM2  [79].  These  compounds  were  designed  on  the 
basis  of  a  spiro-oxindole  core  structure,  which  was  identified 
as  a  scaffold  from  a  substructure  search  of  natural  products. 
Figure  15  shows  the  structures  and  the  competitive  binding 
activities  of  representative  compounds  optimised  on  the 
basis  of  the  spriro-oxindole  template.  The  spiro-oxindoles 
are,  therefore,  exciting  new  additions  of  small  molecules  to 


the  field.  Their  co-crystal  structures  with  MDM2  and 
in  vivo  efficacy  studies  are  eagerly  awaited. 

2.8  a-Helix  mimetic^oKipounds 

A  peptidomimetic-approach  in  general  can  result  in  different 
classes  of  small  molecules.  Protein  binding  interfaces  are  largely 
comprised  of  P-sheets,  a-helices,  turns  and  secondary  loop 
structures.  It  is  possible  to  design  inhibitors  of  protein-protein 
interaction  mimicking  each  of  these  structural  patterns  [80]. 
^■^reviously,  several  compounds  were  successfully  designed  to 
i^^emulate  p-turn  conformations,  which  is  defined  by  the  (j)  and 
\|/  torsion  angles  of  the  t+1  and  z'+2  residues  occupying  the 
turn  region  [80-85] .  Recently,  an  a-helix  mimetic  was  reported 
based  on  terphenyl  derivatives  as  inhibitors  of  the  interaction 
between  calmodulin  (CaM)  and  small  muscle  myosin  light 
chain  kinase  [86].  Such  terphenyl  scaffolds  were  recently 
reported  to  interfere  with  the  p53-MDM2  interaction  [45]. 

Several  novel  inhibitors  of  the  p53— MDM2  interaction 
were  designed  as  a-helix  mimics  using  structure-based  com¬ 
puter-aided  drug  design  methods.  The  CAVEAT  software 
program  was  applied  to  search  the  Available  Chemical 
Directory  and  other  databases  to  identify  lead  compounds 
having  nearly  similar  geometrical  relationship  as  the  Ca-CP 
of  the  triad  motif  residues  (FI 9,  W23,  and  L26)  of  MDM2- 
binding  p53  peptide.  The  identified  hits  were  further  evalu¬ 
ated  by  surface  docking  procedure,  and  the  synthetically 
accessible  scaffolds  were  chosen  for  library  synthesis 
(Figure  16).  Furthermore,  the  procedure  in  general  could  be 
used  to  design  small  molecule  inhibitors  to  interfere  with 
the  protein-protein  interactions  that  are  largely  comprised 
of  helices  of  each  protein  [133].  The  compounds  obtained 
from  this  approach  have  a  variety  of  structural  templates  for 
lead  optimisation. 


3.  Discovery  of  small  molecule  ubiquitin  ligase 
inhibitors 

Another  approach  to  stabilise  the  p53  cellular  level  is  by 
ubiquitin  ligase  inhibitors  that  protect  p53  from  MDM2- 
mediated  ubiquitylation.  The  ubiquitylation  of  target 
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71 


72 


Figure  11.  Representative  substituted  piperizine  inhibitors  of  p53-MDM2. 


proteins  is  regulated  by  the  enzymatic  activity  of  three  pro¬ 
teins.  They  are  the  El  ubiquitin-activating  enzymes,  the  E2 
ubiquitin-conjugating  enzymes  and  E3  ubiquitin-protein 
ligases.  It  has  been  proposed  that  E3  functions  as  ‘docking 
proteins’  for  specific  substrate  protein  and  E2  binding,  and 


that  ubiquitin  is  then  transferred  directly  from  E2s  to  sub¬ 
strates  [87,88].  MDM2  belongs  to  the  second  subclass  of  sin¬ 
gle  subunit  E3  ligase  that  have  the  RING  finger  domain  and 
substrate  recognition  domain  on  a  single  subunit,  whereas 
the  other  subclass  contains  the  RING  finger  domain  and  the 
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Figure  12.  Representative  substituted  piperidine  inhibitors  of  p53-MDM2  interaction. 
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Figure  13.  Representative  aryl  boronic  acid  inhibitors  of  p53-l\/IDI\/l2  interaction. 


Table  5.  Growth  inhibition  of  a  panel  of  human  breast 
cancer  cells  by  compound  88. 


Cell  line 

IC50  (mM) 

MCF7 

1.9 

MDA-MB-231 

1.3 

MCF-10A 

2.5 

MCF-12A 

4.4 

Table  6.  Selective  growth  suppression  of  vl/ild-type  p53- 
expressing  cells  by  fused  indole  inhibitors  of 


p53-MDM2  interaction. 


Compound 


IC50  (mM) 


92 

93 

94 

95 


4^^ 


6p53+ 


HCT116p53-'- 


3.1 
1.8 
<  1 

1.1 
10 


17 

17 

5.6 

43 

12 

48 


substrate  recognition  domain  on  different  subunits  [89].  A 
recent  study  showed  that  small  molecules,  such  as  benzene- 
sulfonamides,  ureas  and  imidazolone  derivatives,  inhibit 
MDM2-mediated  p53  ubiquitylation  (Figure  17)  [90]. 

Another  set  of  compounds,  including  5-deazaflavin,  were 
shown  to  inhibit  MDM2  autoubiquitylation  and  clearly  lead 
to  the  activation  of  p53  function  despite  allowing  the 
stabilisation  of  both  MDM2  and  p53  [9i]. 


5-Deazaflavin  compounds  were  shown  to  effectively  regu¬ 
late  p53  activity  and  anticancer  activities  by  inhibiting 
MDM2-regulated  ubiquitylation  [134].  This  finding  demon¬ 
strated  that  p53  levels  can  be  stabilised  in  mammalian  cells 
after  the  treatment  of  7-nitro-5-deazaflavin.  Three  representa¬ 
tive  compounds  are  shown  in  Figure  18.  Biological  assays  dem¬ 
onstrated  that  the  compounds  only  increased  the  amount  of 
MDM2  and  p53  but  not  p21,  indicating  they  are  specific  for 
E3  activity  of  MDM2. 

A  recent  invention  provided  a  group  of  compounds  as 
inhibitors  of  ubiquitin  ligases  [135].  The  molecules  can  also  be 
used  in  treating  different  diseases  where  ubiquitylation  is 
involved.  Because  MDM2  mediates  multiple  monoubiquityl- 
ation  of  p53  by  a  mechanism  requiring  enzyme  isomerisation, 
such  compounds  can  stabilise  p53  cell  levels.  Several 
representative  compounds  are  shown  in  Figure  19. 

4.  Small  molecule  activators  of  mutant  p53 

The  previously  mentioned  compounds  are  aimed  at  restoring 
or  stabilising  the  wild-type  p53  level  to  control  cell  growth 
and  death.  However,  considering  the  fact  that  over  half  of  all 
cancers  carry  p53  mutations,  the  reactivation  of  p53-specific 
DNA  binding  is  very  important  to  trigger  p53-mediated 
apoptosis  in  cancer  cells  under  pathological  conditions. 
Accordingly,  such  strategies  for  the  treatment  of  mutant  p53- 
carrying  cancers  should  greatly  improve  clinical  outcome  [92]. 
Several  studies  were  earlier  reviewed  to  provide  methods  to 
restore  normal  function  of  mutant  p53  [93].  Some  small  mole¬ 
cules  were  identified  by  a  chemical  library  screening  showing 
anticancer  activity  through  their  ability  to  restore  the  DNA- 
binding  activity  of  a  mutant  p53  (Figure  20)  [94].  Separate 
chemical  library  screening  from  the  National  Cancer  Institute 
discovered  a  novel  compound,  PRJMA-1,  to  induce  mutant 


180 


Expert  Opin.  Then  Patents  (2006)  16(2) 


Deng,  Dayam  &  Neamati 


Figure  15.  Representative  spiro-oxindoles. 

,9 


Figure  16.  Representative  a-helix  mimetic  inhibitors  of  p53-MDM2  interaction. 


Expert  Opin.  Then  Patents  (2006)  16(2) 


181 


Patented  small  molecule  inhibitors  of  p53-MDM2  interaction 


Figure  19.  Arylsulfonamide  and  amide  ubiquitin  ligase  inhibitors. 
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Figure  20.  Examples  of  compounds  reactivating  the  DNA  binding  activity  of  mutant  p53. 
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Figure  21.  Representative  compounds  restoring  apoptosis  inducing  function  of  mutant  p53. 


Table  7.  Calculated  physicochemical  properties  of  some  of  the  patented  inhibitors  of  p53-MDM2  interaction. 


Compound 


Calculated  physicochemical  properties 


MW 


S+logP 


HBA 


HBD 


RB 


3 

8 

9 

12 

18 

31 

42 

52 

54 

59 

61 

65 

66 
6XA’ 


73 

78 

86 

92 

94 

95 
99 


497.41 

639.61 
605.59 
631.56 
413.36 
473.39 

580.29 
600.44<^ 
49?^' 

651.31 

706.44 

360.45 
812.26 
859.06 

532.61 

443.91 

385.92 
535.98 

315.32 

349.30 
516.50 
586.72 


5.49 
5.40 
6  32 

4.05  O' 
6.04 


v5^'6.77 


4.78 

4.57 

4.94 

5.89 

5.6 

2.92 

6.49 

6.32 

3.35 

5.91 

5.10 

5.13 

2.64 

-1.88 

3.69 

8.40 


5 
8 

6 

9 
2 
3 
6 
6 
6 
6 

7 
6 

10 
12 

8 
3 

3 

5 

6 
10 
7 

4 


7 
11 

8 
9 

4 
6 

5 

6 
6 
9 

7 

5 

20 

21 

8 

7 

8 
12 

6 
5 
5 

12 


HBA:  Number  of  hydrogen-bond  acceptors;  HBD:  Number  of  hydrogen-bond  donors;  MW:  Molecular  weight;  RB:  Number  of  rotatable  bonds;  S+logP:  Logarithm  of 
the  octanol-water  partition  coefficient  (calculated  using  ADMET  Predictor  [Simulations  Plus,  Inc.]). 
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p53-dependent  apoptosis  and  inhibit  in  vivo  tumour  growth 
[95].  Furthermore,  PRIMA-1  was  synergistic  with  cisplatin  in 
inducing  apoptosis  [96],  Recently,  maleimide  analogues  were 
reported  to  induce  apoptosis  in  human  tumours  by  restoring 
the  normal  function  of  mutant  p53  [92,136],  The  structures  of 
representative  compounds  are  shown  in  Figure  21. 

5.  Expert  opinion 

p53  is  a  transcription  factor  forming  an  autoregulartory  loop 
with  MDM2.  Over  50%  of  all  cancers  have  mutant  p53,  indi¬ 
cating  its  critical  role  in  the  cell.  Therefore,  restoring  or  stabil¬ 
ising  the  p53  cellular  level  and  its  normal  function  has  become 
one  of  the  major  targets  for  cancer  therapy.  The  disruption  of 
the  p53-MDM2  interaction  has  been  validated  as  a  novel  drug 
target  for  the  development  of  cancer  therapeutics.  In  particu¬ 
lar,  the  discovery  of  the  Nutlin  series  of  compounds  provided 
compelling  evidence  that  the  p53-MDM2  interaction  can  be 
successfully  manipulated  with  small  molecule  inhibitors.  The 
success  of  these  compounds  also  provides  strong  support  for 
the  therapeutic  application  of  selective  disruption  of  protein- 
protein  interactions.  Protein— protein  interactions  consist  of 
large  surfaces  and  multiple  interactions  of  diverse  nature. 
Many  of  the  inhibitors  of  the  p53— MDM2  interaction 
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discussed  in  this  review  contain  several  hydrophobic  groups 
and  have  a  high  molecular  weight.  To  put  these  compounds 
into  perspective  the  authors  have  calculated  drug-like  proper¬ 
ties  (Lipinski’s  rule-of-five)  of  some  of  the  selected  compounds 
from  each  class  of  patented  inhibitors  (Table?).  Interestingly, 
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functional  groups.  These  inhibitors  represent  a  novel  cla^g"^. 
therapeutics  targeted  at  protein-protein  interactions:  upon 
their  in  vivo  success  some  modifications  to  Lipinski’s  rule-of- 
five,  which  was  formulated  on  the  basis  of  properties  of  drugs 
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The  human  epiderm^  growth  factor  receptor  (HER)  family  is  a  highly 
explored  and  promismg  anticancer  drug  target.  At  present,  several  investi¬ 
gational  agents  targeted  to  the  HER  family  of  receptors  are  in  various 
stages  of  development.  Five  drugs  are  already  in  the  clinic  for  the  treat¬ 
ment  of  cancers  that  overexpress  HER  family  receptors.  Two  FDA-approved 
small-molecule  drugs,  gefitinib  and  erlotinib,  inhibit  HER1  tyrosine  kinase 
activity.  Two  mAbs,  cetuximab  and  panitumumab,  target  the  extracellular 
domain  of  HER1,  and  another,  trastuzumab,  targets  the  extracellular 
domain  of  HER2.  HER2  is  a  prominent  member  of  the  HER  family  of  recep¬ 
tor  tyrosine  kinases  and  serves  as  a  preferred  dimerization  partner  for 
other  HER  family  members.  This  paper  reviews  recently  patented 
small-molecule  inhibitors  of  HER2  receptor  kinase  activity,  and  inhibitors  of 
HER2  expression  and  shedding.  Apart  from  the  well-explored  quinazoline 
class  of  compounds  (e.g.,  lapatinib),  arylazole,  benzodithiazole,  pyrrolopy- 
ridazine,  pyrrolotriazine  and  pyrrolopyrimidine  classes  of  compounds  were 
also  claimed  as  HER2  tyrosine  kinase  inhibitors.  Most  of  these  compounds 
show  considerable  activity  against  all  the  HER  family  as  well  as  members 
from  different  families  of  tyrosine  kinases.  It  remains  to  be  established 
how  the  combination  of  selective  HER  inhibitors  compare  with  the  sin¬ 
gle-agent  pan-kinase  inhibitors  in  disrupting  HER  family  mediated  signal¬ 
ling  pathways.  Such  information  is  of  paramount  importance  in  the  clinical 
development  of  HER-targeted  inhibitors. 


Keywords:  anticancer  drug,  HER2,  monoclonal  antibody,  small-molecule, 
tyrosine  kinase  inhibitor 

Expert  Opin.  Ther.  Patents  (2007)  17(1):83'102 

1.  Introduction 


Molecularly  targeted  anticancer  drugs  represent  a  rational  advancement  in  the  mod¬ 
ern  drug  discovery  paradigm.  In  addition  to  cytotoxic  and  DNA-targeted  anticancer 
agents,  a  large  number  of  mechanism-based  molecular  targets  are  under  intensive 
investigation.  Clinical  success  of  several  molecularly  targeted  therapeutics  has  accel¬ 
erated  tremendous  interest  in  this  field.  The  human  epidermal  growth  factor  recep¬ 
tor  (HER,  ErbB)  family  is  an  excellent  example  for  mechanism-based  anticancer 
molecular  targets.  The  HER  family  consists  of  four  closely  related  transmembrane 
proteins,  epidermal  growth  factor  receptor  (EGFR/ErbBl/HERl),  human 
epidermal  growth  factor  receptor-2  (HER2/Neu/ErbB2),  HER3/ErbB3  and 
HER4/ErbB4.  The  HER  family  receptors  are  among  the  most  studied  cell  signalling 
families  in  cancer  biology  [i].  Each  protein  has  an  extracellular  ligand-binding 
domain,  a  transmembrane  region  and  an  intracellular  cytoplasmic  domain  [2].  The 
extracellular  domain  is  involved  in  recognizing  and  binding  ligands  that  are  able  to 
activate  receptor.  The  cytoplasmic  domain  contains  a  tyrosine  kinase  region  and  a 
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carboxy  terminal  tail  that  harbors  tyrosine  autophosphoryla¬ 
tion  sites.  The  HER  family  receptors  have  similar  structures 
and  a  high  degree  of  homology  in  the  tyrosine  kinase  region, 
although  they  possess  diverse  characteristics  that  determine 
their  signalling  specificity.  The  tyrosine  kinase  domains  of 
HER2  and  -4  show  ~  80%  homology  to  that  of  HERl, 
whereas  HER3  lacks  intrinsic  tyrosine  kinase  activity  [3].  Inter¬ 
estingly,  the  extracellular  domains  are  less  conserved  among 
the  four  receptors,  which  is  indicative  of  different  specificity  in 
ligand  binding. 

The  HER  receptors  exist  as  inactive  monomers.  A  variety 
of  epidermal  growth  factor  ligands  can  bind  to  extracellular 
domains  of  HERl,  -3  and  -4  and  activate  these  receptors. 
These  ligands  have  different  specificities  for  each  receptor, 
resulting  in  different  cellular  signalling  effects  [4,5].  There  is  no 
known  ligand  for  HER2,  rendering  it  an  orphan  receptor.  On 
binding  to  ligands,  the  extracellular  domains  of  HERl,  -3  and 
-4  receptors  undergo  conformational  changes  that  allow  acti¬ 
vation  of  the  receptors  by  either  a  homodimerization  or  a  het- 
erdimerization  process,  which  is  essential  to  stimulate  the 
tyrosine  kinase  activity  of  the  receptors  for  subsequent  genera¬ 
tion  of  an  intracellular  signal.  The  pairing  of  two  receptors  of 
the  same  type  results  in  homodimers,  whereas  pairing  of  two 
different  type  of  receptors  produces  heterodimers.  Interest¬ 
ingly,  constitutive  conformation  of  the  extracellular  domain 
of  HER2,  an  orphan  receptor,  resembles  the  conformation  of 
ligand  activated  extracellular  domains  of  the  other  HER 
receptors  [6,7].  This  autoactivated  conformation  enables 
HER2  to  perform  as  a  preferred  dimerization  partner  for  all 
the  other  HER  receptors.  It  has  been  observed  that  the  rela¬ 
tive  expression  levels  of  HER  receptors  and  the  ligands  dictate 
a  receptors’  ability  to  form  either  a  ‘homo’  or  a  ‘hetero’  dimer. 
It  also  has  been  demonstrated  that  HER  receptors  can  acquire 
different  signalling  properties  following  dimerization  with  dif¬ 
ferent  partners.  For  example,  a  heterodimeric  complex  con¬ 
taining  HER2,  and  an  inactive  tyrosine  kinase  receptor, 
HER3,  can  form  the  most  potent  HER  signalling  complex  in 
terms  of  cell  proliferation  and  transformation  [8]. 

HER  receptor  dimerization  induces  tyrosine  kinase  cata¬ 
lytic  activity,  which  leads  to  autophosphorylation  of  several 
key  tyrosine  residues  in  the  activation  loop  of  the  tyrosine 
kinase  domain  within  the  receptor’s  carboxyl-terminal  tail  via 
transfer  of  y-phosphates  from  bound  ATP.  The  resulting  phos- 
photyrosine  residues  act  as  docking  sites  for  a  number  of  sig¬ 
nal-transducing  enzymes  and  adaptor  proteins  [i].  At  least  two 
major  pathways  are  involved  in  HER  signalling,  the  MAPK 
and  phosphatidylinositol  3-kinase  (PI3K)/Akt  pathways 
(Figure  1)  [9,10].  The  four  HER  receptors  share  most  of  the 
same  intracellular  second  messengers  generated  when  they  are 
activated  [ii].  These  include  adaptor  proteins  such  as  She  and 
Grb2,  which  couple  the  HER  receptors  to  the  MAPK  path¬ 
way,  kinases  such  as  c-src  and  PI3K,  and  protein  tyrosine 
phosphatases  such  as  SHPl  and  -2  [ii].  The  various  postrecep¬ 
tor  downstream  signalling  cascades  activated  upon  ligand 
binding  and  receptor  dimerization  result  in  different  cellular 


effects.  For  example,  the  MAPK  pathway  is  important  for  cell 
proliferation,  whereas  the  PI3K  pathway  has  a  predominant 
role  in  cell  survival  [12,13].  The  HER  receptors  and  their  lig¬ 
ands  are  frequently  amplified  in  various  human  cancers  and 
play  an  important  role  in  the  pathogenesis  of  these  diseases. 

1.1  The  HER  family:  a  rational  anticancer  drug  target 

In  a  variety  of  human  cancer  cells,  aberrant  signalling  involving 
HER  receptors  stimulates  pathways  that  activate  many  of  the 
properties  associated  with  cancer,  including  proliferation,  migra¬ 
tion,  metastases,  angiogenesis  and  resistance  to  apoptosis.  Due  to 
the  higjjLpreqiiency  of  abnormalities  in  receptor  signalling  in 
humaaiancers,  the  HER  family  is  an  attractive  target  for  thera¬ 
peutic  development.  The  frequently  observed  causes  of  signalling 
abnormalities  in  cancers  involving  HER  receptors  are  overexpres¬ 
sion  of  receptors,  overproduction  of  growth  factor  ligands  and 
ligand-independent  receptor  activation.  The  third  condition  is 
related  to  a  mutant  receptor,  EGFRvIII,  which  is  commonly 
detected  in  many  types  of  cancers.  EGFRvIII  lacks  a  major  seg¬ 
ment  in  the  extracellular  ligand  binding  domain.  This  mutation 
results  in  ligand-independent,  constitutive  activation  of  the 
receptor  [i4].  The  HER  receptors  are  overexpressed  or  dysregu- 
lated  in  a  wide  variety  of  cancers  including  breast,  colorectal, 
ovarian,  prostate  and  non-small-ceU  lung  cancers.  Importantly, 
overexpression  of  HER  receptors  is  associated  with  poor  disease 
prognosis  and  reduced  survival.  Mounting  preclinical  and  clini¬ 
cal  evidence  supports  the  rationale  behind  the  HER  family  tar¬ 
geted  anticancer  therapeutic  approaches.  Although  the 
complexities  of  the  signalling  pathways  involving  the  HER  fam¬ 
ily  of  receptors  are  not  fully  understood,  several  possible  points 
within  the  pathways  have  been  identified  to  be  therapeutically 
advantageous  if  interrupted  [15].  The  most  promising  and 
advanced  therapeutic  strategies  targeting  HER  receptors  are 
mAbs  and  small-molecule  tyrosine  kinase  inhibitors.  Intensive 
research  over  the  last  20  years  in  these  areas  has  yielded  a  number 
of  promising  anticancer  therapeutics.  As  of  today,  five  anticancer 
therapeutics  targeting  the  HER  receptor  family  have  been 
approved  by  the  US  FDA.  Two  of  them  are  small-molecule  tyro¬ 
sine  kinase  inhibitors,  gefitinib  (Figure  2)  and  erlotinib  (Figure  2), 
targeted  to  the  HERl  tyrosine  kinase  domain.  The  remaining 
three  are  mAbs.  Trastuzumab,  a  humanized  mAh,  is  targeted  to 
the  extracellular  domain  of  HER2.  Cetuximab  is  also  a  human¬ 
ized  mAh  targeted  to  the  extracellular  domain  of  HERl. 
Recently,  the  US  FDA  approved  panitumumab,  the  first  fully 
human  mAb  for  the  treatment  of  patients  with  HERl -expressing 
metastatic  colorectal  cancer.  Like  cetuximab,  panitumumab  is 
also  targeted  to  the  extracellular  domain  of  HERl .  At  present, 
several  small-molecule  HER  tyrosine  kinase  inhibitors  and  HER 
extracellular  domain-targeted  antibodies  are  being  investigated  in 
various  stages  of  clinical  trials.  Considering  the  enormous 
research  activity  in  the  discovery  of  anticancer  therapeutics  tar¬ 
geting  the  HER  receptor  family,  the  emphasis  of  this  review  is 
restricted  to  HER2-targeted  small-molecule  therapeutics. 
Recently  patented  small-molecule  inhibitors  of  tyrosine  kinase 
activity  and  expression  and  shedding  of  HER2  are  reviewed. 
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Figure  1.  Schematic  representation  of  the  HER  family  mediated  signaling  pathways.  HER3  lacks  Intracellular  tyrosine  kinase 
activity;  however,  it  acts  as  a  dimerization  partner  for  other  members  of  the  HER  family  of  receptors. 
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Figure  2.  The  US  FDA  approved  small-molecule  FIER  tyrosine  kinase  inhibitors. 


2.  Anticancer  therapeutics  targeting  HER2 

HER2  is  a  unique  member  of  the  HER  receptor  family.  As 
discussed  above,  HER2  is  an  orphan  receptor  and  operates  as 
a  preferred  heterodimeric  partner  to  the  other  HER  recep¬ 
tors.  The  HER2  extracellular  domain  exists  in  a  heterodimer- 
ization-favorable  conformation  on  the  cell  surface,  which 
underlies  its  enhanced  capacity  for  heterodimerization  with 
other  HER  family  members.  Although  both  homo-  and  het- 
rodimerization  stimulate  HER  signalling  cascades,  het¬ 
erodimers  are  more  potent  and  mitogenic.  HER2  is 
overexpressed  in  25  -  30%  of  breast  cancers  and  high  expres¬ 
sion  levels  correlate  with  poor  prognosis  in  cancers  of  breast, 
ovary,  colon  and  bladder  [16-18].  HER2  has  become  an  impor¬ 
tant  anticancer  therapeutic  target,  particularly  for 


o 

other  HER  family  members.  A  recently  completed  Herceptin 
Adjuvant  (HERA)  trial  study  concludes  that  treatment  with 
trastuzumab  aher  adjuvant  chemotherapy  significantly  improves 
disease-free  survival  among  women  with  HER2-positive  breast 
cancer  [23].  To  maximize  trastuzumab-based  treatment  in 
patients  with  HER-positive  cancers,  novel  combinations  that 
include  tyrosine  kinase  inhibitors  and  antibodies  targeting  other 
HER  receptors  should  be  explored. 

Pertuzumab  is  a  second  antibody  targeted  to  the  HER2 
extracellular  domain.  It  disrupts  HER2-mediated  signalling 
by  sterically  blocking  a  binding  pocket  in  the  extracellular 
domain,  which  is  required  for  HER2  heterodimerization 
with  other  members  of  the  HER  family  [24,25].  Several  pre- 
clinical  studies  show  that  pertuzumab  is  as  active  as  trastuzu¬ 
mab  on  breast  cancer  cells  expressing  high  levels  of  HER2. 
However,  unlike  trastuzumab,  it  is  also  effective  in  inhibiting 


HER2-expressing  primary  and  metastatic  breast  cancers.  The 
first  clinically  used  HER2  specific  anticancer  therapeutic  is  ^the  growth  of  cancer  cells  that  do  not  overexpress  the  HER2 


trastuzumab,  a  recombinant  humanized  mAb  targeted  at  the 
extracellular  domain  of  HER2.  Several  antibodies  and 
small-molecule  tyrosine  kinase  inhibitors  targeting  HER2  are 
in  various  stages  of  development  at  present. 

2.1  Monoclonal  antibodies  targeting  the  HER2  C 
extracellular  domain 


receptor.  Its  efficacy  in  cancers  with  low  HER2  expression 
levels  could  be  attributed  to  its  ability  to  disrupt  HER2  het¬ 
erodimerization  with  other  HER  receptors.  Pertuzumab  is 
presently  in  Phase  II  clinical  evaluation  for  its  effect  on 
patients  with  HER2-overexpressing  advanced  or  metastatic 
breast  cancers  [26[. 


mAbs  are  highly  specific  therapeutics.  Several  antibodies  tar¬ 
geted  at  HER  family  members  are  in  development  for  treat¬ 
ment  of  various  cancers.  At  present,  trastuzumab  is  the  only 
therapy  approved  for  the  treatment  of  HER2-positive  cancers. 
A  recent  X-ray  crystal  structure  reveals  several  critical  interac¬ 
tions  between  the  HER2  extracellular  domain  and  trastuzu¬ 
mab.  Trastuzumab  shows  both  cytostatic  and  cytotoxic 
activities  [19].  It  inhibits  tumor  growth  by  accelerating '  the 
internalization  and  degradation  of  the  HER2  receptor  [20,21]. 
It  has  been  also  shown  that  trastuzumab  has  antiangiogenic 
activity  and  downregulates  proteins  involved  in  angiogenesis 
[22].  It  also  has  been  demonstrated  that  trastuzumab  shows 
synergism  in  combination  with  a  variety  of  chemotherapeutic 
agents  in  HER2  positive  cancers.  However,  the  majority  of 
cancers  that  initially  respond  to  trastuzumab  treatment  start 
to  progress  over  the  course  of  a  year.  Some  of  the  reasons  for 
resistance  to  trastuzumab  include  altered  HER2-trastuzumab 
interactions  and  increased  compensatory  cell  signalling  by 


2.2  Small  molecules  targeting  the  HER2  tyrosine 
kinase  domain 

Small-molecule  tyrosine  kinase  inhibitors  targeted  to  the 
HER  receptor  family  act  at  an  intracellular  level.  These  com¬ 
pounds  compete  with  ATP  for  binding  to  the  ATP  site  on  the 
receptor’s  intracellular  tyrosine  kinase  domain.  By  binding  to 
the  ATP  site  and  inhibiting  catalytic  activity,  they  block  auto- 
pbosphorylation  of  key  tyrosine  residues  on  the  receptor 
kinase  domain.  This  disrupts  activation  of  the  receptor  medi¬ 
ated  downstream  signalling  cascades.  Several  small-molecule 
HER  tyrosine  kinase  inhibitors  are  in  clinical  development. 
At  present,  there  are  no  FDA-approved  tyrosine  kinase  inhibi¬ 
tors  selectively  targeting  HER2  available  for  cancer  treatment. 
However,  the  two  clinically  approved  HER  1 -specific  tyrosine 
kinase  inhibitors  gefitinib  and  erlotinib  also  inhibit  HER2 
phosphorylation  and  downstream  signalling  cascades  in 
HER2/HER3-overexpressing  cancer  cells  [27-29].  It  has  been 
observed  that  most  of  the  HER  tyrosine  kinase  inhibitors  in 


86 


Expert  Opin.  Then  Patents  (2007)  17(1) 


Dayam,  Grande,  Al-Mawsawi  &  Neamati 


Table  1.  Growth  inhibitory  activities  of  representative  arylalkoxy  oxazole  and  thiazole-based  HER  tyrosine  kinase 
inhibitors  in  A549  cells.  A 
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Table  2.  Inhibition  of  HER2  phosphorylation  by  representative  pentaflurosulfanyl  derivatives  of 
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97.5 


91.9 
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clinical  development  also  show  cross  receptor  activity.  This 
phenomenon  could  be  attributed  to  the  high  homology 
present  in  kinase  domains  of  HER  receptors.  At  this  point  of 
time,  it  is  not  clear  that  HER  tyrosine  kinase  inhibitors  that 
show  either  a  broad  activity  spectrum  (those  that  inhibit  more 
than  one  receptor  simultaneously)  or  receptor  specificity  yield 
an  optimum  therapeutic  response  with  minimum  or  no  toxic 
effects  in  clinical  settings.  Encouraging  results  show  that  clini¬ 
cally  approved  small-molecule  tyrosine  kinase  inhibitors  have 
a  favorable  safety  profile  and  can  be  combined  with  other 
forms  of  chemotherapy  or  radiation  therapy  [30].  Lapat- 
inib/GW-20l6  (Figure  3)  is  a  reversible  dual  HER2  and  -1 
tyrosine  kinase  inhibitor  [31-33].  Lapatinib  also  inhibits 
p95ErbB2  (truncated  HER2  receptor)  kinase  activity  in  vitro 


and  in  tumor  xenografts  [34].  The  combination  of  lapatinib 
with  trastuzumab  showed  synergistic  effects  on  HER2- 
expressing  and  trastuzumab-treated  cancer  cell  lines  [35].  The 
combination  of  lapatinib  with  trastuzumab  enhances  apopto¬ 
sis  in  HER2-overexpressing  breast  cancer  cell  lines  [36].  A 
recently  completed  Phase  III  clinical  trial  showed  that  lapat¬ 
inib  significantly  improves  disease  progression-free  survival  in 
patients  with  HER2-positive  advanced  and  metastatic  breast 
cancers  whose  disease  had  progressed  following  treatment 
with  trastuzumab  and  other  cancer  therapies  [37].  Regulatory 
approval  from  the  FDA  for  clinical  use  of  lapatinib  is  awaited. 
This  quinazoline  derivative  represents  a  first-generation 
HER2-targeted  small-molecule  anticancer  agent.  Caner- 
tinib/CI-1033  (Figure  3)  is  an  irreversible  pan-HER  tyrosine 
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4  CI-1033 


Figure  3.  Structure  of  lapatinib  and  CI-1033,  pan-HER  tyrosine  kinase  inhibitors  currently  in  clinical  trials. 
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IC50  concentrations 
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range 


of  7  -  70  nM. 


kinase  inhibitor  [38-41].  CI-1033  inhibits  HERl,  -2  and  -4 
kinase  activities  with  similar  potency.  Phase  II  clinical  trials  of  Compound  16  containing  an  arylalkoxy  oxazole  core  inhibited 

CI-1033  are  under  way.  Recently  a  number  of  small  mole-  cell  growth  at  a  concentration  of  7  nM. 

cules  with  diverse  structural  scaffolds  have  been  claimed  as  Several  pentaflurosulfanyl  derivatives  of  arylalkoxy  oxazole 
HER2  tyrosine  kinase  inhibitors  in  several  patents.  Recent  core-containing  compounds  show  significant  inhibition  of 

advances  in  the  design  and  discovery  of  small-molecule  HER2Q  HER2  phosphorylation  at  a  tested  concentration  of  1  pM 

tyrosine  kinase  inhibitors  are  reviewed  below. 


-u 

2.2.1  Quinazolines  as  HER  tyrosine  kinase  inhibitors 
Quinazolines  have  been  one  of  the  most  highly  exploraST 
chemical  classes  in  the  discovery  of  HER  tyrosine  kiif»& 
inhibitors.  Clinically  used  HERl  tyrosine  kinase  inhibitors 
gefitinib  and  erlotinib,  and  the  most  advanced  HER2-tar- 
geted  inhibitors,  lapatinib  and  CI-1033,  belong  to  the 
quinazoline  class  of  compounds.  These  compounds  have 
been  designed  to  compete  with  ATP  by  binding  at  the  ATP 
binding  site  on  the  HER  tyrosine  kinase  domain.  Several 
novel  quinazoline  HER  tyrosine  kinase  inhibitors  were  dis¬ 
closed  in  a  recent  patent  [101].  Structures  of  some  representa¬ 
tive  quinazolines  are  shown  in  Figure  4.  Although  the 
quinazoline  core  has  been  highly  explored  in  the  design  of 
HER  tyrosine  kinase  inhibitors,  these  novel  compounds  pos¬ 
sess  a  unique  hydrazinecarboximidamide  substitution  at  the 
sixth  position  of  the  quinazoline  core,  thus  making  them 
different  from  previously  reported  quinazolines.  These  com¬ 
pounds  inhibit  HERl,  -2  and  -4  tyrosine  kinase  inhibitors 
[101].  They  are  also  being  claimed  as  serine,  threonine  and 
dual  specific  kinase  inhibitors.  However,  details  of  receptor 
specific  activities  and  in  vivo  anticancer  activity  profiles  of 
compounds  related  to  this  patent  are  not  described.  Com¬ 
pounds  presented  in  this  patent  are  reported  to  inhibit  HER 


tyrosine  kinase  activity  with 
range  of  <  1  nM  —  50  pM. 
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2.2.2  Arylazole-based  HER2  tyrosine  kinase  inhibitors 
Compounds  containing  an  arylalkoxy  oxazole  or  thiazole  core 
are  potent  inhibitors  of  HER2  tyrosine  kinase  activity  [102]. 
Structures  of  selected  compounds  along  with  their  antiprolifer¬ 
ative  activities  in  A549  (human  lung  carcinoma)  cells  are  given 
in  Table  1.  These  compounds  inhibited  A549  cell  proliferation 


[103].  The  per  cent  inhibition  of  these  compounds  on  HER2 
phosphorylation  is  shown  in  Table  2.  Arylazole  (27)  inhibited 
97.5%  of  HER2  phosphorylation  at  a  concentration  of  1  pM, 
and  is  the  most  potent  HER2  tyrosine  kinase  inhibitor  from 
this  series  of  compounds.  Several  structurally  related  com¬ 
pounds  with  the  arylalkoxy  oxazole  core  were  also  claimed  as 
HER2  tyrosine  kinase  inhibitors  [104].  These  compounds 
show  significant  inhibition  of  HER2  phosphorylation  in 
Calu3  tumor  cells.  Some  of  the  selected  compounds  inhibited 
70  -  86%  of  HER2  phosphorylation  at  1  pM.  Structures  of 
representative  compounds  and  their  HER2  tyrosine  kinase 
inhibitory  activities  are  presented  in  Table  3. 

Arylalkoxy  oxazoles,  in  which  a  diazine  moiety  represents 
the  aryl  group,  are  also  potent  inhibitors  of  HER2  tyrosine 
kinase  activity  [105].  Structures  of  three  representative  com¬ 
pounds  and  their  HER2  tyrosine  kinase  inhibitory  activities 
are  presented  in  Table  4.  In  addition,  indole  derivatives  of 
arylazoles  also  inhibit  HER  tyrosine  kinase  and  exhibit  signif¬ 
icant  growth  inhibition  of  HEK293  (human  embryonic  kid¬ 
ney)  cells  [106].  Structures  of  representative  compounds  and 
their  anticancer  activities  are  presented  in  Table  5.  Further¬ 
more,  amide  derivatives  of  arylazoles  also  inhibited  HER2 
phosphorylation  at  1  pM  and  blocked  HEK293  cell  growth 
in  the  nanomolar  range  [107].  Structures  of  selected  derivatives 
and  their  per  cent  inhibition  of  HER2  phosphorylation  are 
presented  in  Table  6.  Previously,  several  structurally  similar 
arylalkoxy  oxazoles  were  reported  to  show  antidiabetic  effects 
[42].  Some  of  these  compounds  show  strong  glucose  lowering 
and  lipid  lowering  effects  in  diabetic  animal  models.  The  anti¬ 
diabetic  effects  of  these  compounds  are  shown  to  be  con¬ 
nected  to  their  potent  agonistic  activity  on  peroxisome 
proliferators-activated  receptor-y  (PPARy).  Several  structur¬ 
ally  related  arylalkoxy  oxazoles  containing  a  2,4-thiazolidine- 
dione  or  2,4-oxazolidinedione  ring  substitution  on  the 
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Figure  4.  Representative  quinazoline-based  HER  tyrosine  kinase  inhibitors. 
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Figure  6.  Representative  benzodithiazole-based  HER  tyrosine  kinase  inhibitors. 
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Table  3.  Inhibition  of  HER2  phosphorylation  by  representative  arylazole-based  HER2  tyrosine  kinase  inhibitors. 


Compound 


Structure 


Percent  inhibition  of  HER2 
phosphorylation  at  1  pM 
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arylalkoxy  group  also  show  antidiabetic  effects  through  their 
potent  agonistic  activity  on  PPARy  [43]. 

TAK-165  (48),  an  arylalkoxy  oxazole-containing  HER2 
tyrosine  kinase  inhibitor,  is  presently  under  clinical  studies 
(Figure  5).  A  Phase  I  study  using  a  once-daily  oral  dosing  regi¬ 
men  of  TAK-165  in  patients  with  advanced  cancers  that 
express  HER2  has  been  completed  [44],  TAK-165  selectively 
inhibits  HER2  tyrosine  kinase  activity  at  an  ICjg  value  of 
6  nM  and  does  not  inhibit  other  types  of  tyrosine  kinases  up 
to  25,000  nM  [45].  It  was  also  shown  to  inhibit  HER2  phos¬ 
phorylation  and  its  downstream  PI3K  and  MAPK  pathways 
in  vitro  in  HER2-overexpressing  BT474  breast  cancer  cells. 
TAK-165  also  inhibited  growth  of  HER2-expressing  bladder, 
kidney  and  androgen-independent  prostate  cancers.  This 
indicates  a  new  hope  for  the  possible  therapeutic  application 
of  HER2  tyrosine  kinase  inhibitors  in  the  treatment  of  can¬ 
cers  other  than  HER2-overexpressing  breast  cancers  and,  in 
particular,  hormone-refractory  prostate  cancer.  Details  of  the 
preclinical  and  clinical  activity  of  the  arylazole  class  of  com¬ 
pounds  are  eagerly  awaited.  It  is  also  interesting  to  note  that 
TAK-165  inhibited  HER2  phosphorylation  in  human  rheu¬ 
matoid  synovial  cells,  indicating  its  therapeutic  applications 
beyond  cancer  treatment  [108[. 

2.23  Benzo-  and  azabenzodithiazole-based  HER 
tyrosine  kinase  inhibitors  4^ 

Several  benzo-  and  azabenzodithiazoles  were  shown  to  inhibit 
HERl,  -2  and  -4  kinases  with  ICjg  values  of  0.001  -  25  pM 
[109].  However,  details  on  receptor  specific  activity  and  anti¬ 
cancer  properties  of  the  claimed  compounds  are  not  disclosed. 
Structures  of  some  representative  compounds  are  shown  in 
Figure  6.  Previously,  these  compounds  were  not  reported  to 


show  significant  biological  activity.  However,  compounds 
containing  benzodithiazole  were  claimed  as  herbicides  and 
plant  growth  regulators  [iio[. 

2.2.4  Pyrroiopyridazine  and  pyrrolotriazine-based  HER 
tyrosine  kinase  inhibitors 

Pyrroiopyridazine  and  pyrrolotriazine  analogs  are  growing 
examples  of  HER  tyrosine  kinase  inhibitors  [111-117].  Selected 
compounds  representing  structural  diversity  are  shown  in 
Figure  7  [111,112].  Pyrrolopyridazines  are  reported  to  inhibit 
HER  tyrosine  kinases  with  IC50  values  of  0.001  —  25  pM. 
Some  of  these  compounds  are  also  shown  to  inhibit  kinase 
activities  of  VEGF  receptor  (VEGFR)-2  and  fibroblast 
growth  factor  receptor  (FGFR)-l  receptors  with  IC50  values  of 
<  80  pM.  Although  several  hundreds  of  compounds  were 
claimed  as  HER  tyrosine  kinase  inhibitors,  details  on  receptor 
selective  inhibitory  activities  and  anticancer  properties  are  not 
described.  Apart  from  blocking  phosphorylation  activity  of 
tyrosine  kinase  receptors,  pyrrolopyridazines  also  inhibit 
activity  of  a  variety  of  enzymes  implicated  in  different  disease 
conditions.  For  example,  pyrrolopyridazine-containing  com¬ 
pounds  were  reported  to  inhibit  the  activity  of  15 -lipoxygen¬ 
ase  (15-LO)  enzyme  [46].  Gonsidering  the  role  of  15-LO  in 
atherosclerosis,  these  compounds  have  potential  therapeutic 
benefit  in  the  treatment  of  atherosclerosis  and  other 
1 5-LO-mediated  disease  conditions.  Several  pyrrolopyri¬ 
dazines  are  claimed  as  corticotropin-releasing  factor  (GRF) 
receptor  antagonists  [118,119].  These  compounds  are  reported 
to  possess  therapeutic  applications  against  CNS  disorders 
such  as  anxiety  and  depression  associated  with  deregulation  or 
hypersecretion  of  the  GRF  receptor.  Pyrroiopyridazine  ana¬ 
logs  also  inhibit  phosphodiesterase-4  (PDE-IV)  and  produc- 
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Table  4.  Inhibition  of  HER2  phosphorylation  by  representative  diazinealkoxy  oxazole-based  HER2  tyrosine  kinase 
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Table  5.  Growth  inhibitory  activities  of  representative  indole  derivatives  of  arylazole-based  HER2  tyrosine  kinase 
inhibitors  in  HEK293  cells. 
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Table  6.  Inhibition  of  HER2  phosphorylation  by  representative  amide  derivatives  of  arylazole-based  HER2  tyrosine 
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tion  of  TNF-a,  which  is  important  for  the  treatment  of 
cancer,  inflammatory  and  autoimmune  diseases  [120]. 
Recently,  several  pyrrolopyridazine  core  containing  com¬ 
pounds  were  reported  as  MAPK/ERK  (extracellular 
signal-regulated  kinase)  (MEK)  inhibitors  [471.  Some  of  the 
selected  compounds  inhibited  MEK  activity  at  <  1  pM  and 
show  cell  growth  inhibitory  activity  on  human  tumor  cell 
lines  HT29  and  Colo205  at  <  8  pM.  In  summary,  the  pyrrol¬ 
opyridazine  core  is  an  excellent  scaffold  in  designing  of  a  vari¬ 
ety  of  drugs,  including  inhibitors  of  EIER  family  receptor 
tyrosine  kinases. 

The  4-aminopyrrolotriazine-containing  compounds  are 
interesting  examples  of  HER  family  receptor  tyrosine  kinase 
inhibitors  (Figure  8)  [113].  Unfortunately,  their  HER  tyrosine 
kinase  inhibitory  activities  and  anticancer  properties  have  not 
been  reported.  The  4-aminopyrrolotriazine  core  has  been 


considered  as  a  novel  template  for  the  design  of  receptor  tyro¬ 
sine  kinase  inhibitors  [48].  Several  compounds  containing  the 
4-aminopyrrolotriazine  core  have  been  reported  as  dual  inhib¬ 
itors  of  HERl  and  -2  kinase  inhibitors  [49[.  For  example, 
BMS-599626  (72)  is  a  potent  inhibitor  of  HERl  and  -2 
kinase  activity  [50].  BMS-599626  also  inhibited  proliferation 
of  a  collection  of  diverse  cancer  cell  lines  that  are  dependent 
on  signalling  from  either  HERl  or  -2.  However,  how 
BMS-599626  inhibits  HERl  and  -2  kinase  activity  and  sig¬ 
nalling,  and  the  mechanism  of  binding  to  the  receptors  is  not 
clearly  understood.  BMS-599626  acted  as  an  ATP  com¬ 
petitive  inhibitor  in  HERl  inhibition.  Interestingly, 
BMS-599626  inhibited  HER2  with  an  ATP  noncompetitive 
mechanism.  Recently  reported  preclinical  results  project 
BMS-599626  as  a  promising  candidate  for  further  clinical 
development  for  cancers  with  HERl  and  -2  overexpression. 
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Figure  7.  Representative  pyrrolopyridazine-based  HER  tyrosine  kifl^i^  inhibitors  (continued) 
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Figure  8.  Representative  morpholinylmethyl  ester  derivatives  of  pyrrolotriazine-based  HER  tyrosine  kinase  inhibitors. 
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Figure  9.  Substituted  pyrrolotriazine-base  R  tyrosine  kinase  inhibitors. 
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Selected  structurally  diverse  5-substituted  4-aminopyrrolotri- 
azine  HER  tyrosine  kinase  inhibitors  are  shown  in  Figure  9 
[114],  Some  of  these  compounds  were  reported  to  inhibit 
HERl,  -2  and  -4  kinases  with  ICjq  values  of  0.001  —  25  pM. 
Aminocyclohexyloxy  derivatives  of  4-aminopyrrolotriazines 
and  other  analogs  are  claimed  as  HER  tyrosine  kinase  inhibi¬ 
tors  [115,116],  Structures  of  representative  aminocyclohexyloxy 
derivatives  are  shown  in  Figure  10.  These  compounds  also 
inhibit  receptor  tyrosine  kinases  such  as  VEGFR2,  FGFRl 
and  platelet-derived  growth  factor  receptor  in  addition  to 
HER  family  receptor  kinases  [ii6]. 

Several  5-  and  6-substituted  4-aminopyrrolotriazine  core- 
containing  compounds  are  reported  to  inhibit  VEGFR2  and 
FGFRl  kinases  at  <  100  nM  [5i,52].  Selected  compounds  show 
metabolic  stability  and  acceptable  pharmacokinetic  profiles 
after  oral  administration  in  mice  xenografts.  Furthermore, 
some  of  these  compounds  are  proposed  to  bind  at  the  ATP 
binding  site  of  VEGFR2  and  form  H-bonding  interactions 
with  amino  acid  residues  from  the  hinge  region.  Several  ami- 
nopiperidine  and  diazepan  derivatives  of  4-aminopyrrolotri- 
azines  are  also  reported  as  HER  tyrosine  kinase  inhibitors. 
Three  representative  compounds  are  shown  in  Figure  11.  As 
described  above,  4-aminopyrrolotriazine  analogs  inhibit  HER 
tyrosine  kinases  with  ICjq  values  of  0.001  —  25  pM  [117]. 
Moreover,  the  4-aminopyrrolotriazine  core  is  also  considered 
as  a  lead  template  in  the  design  of  p38  kinase  inhibitors  for 
the  treatment  of  inflammatory  disorders  and  potentially  other 
p38  kinase-mediated  diseases  [121]. 

2.2.5  Pyrrolopyrimidine-based  HER2  tyrosine  kinase  ^ 
inhibitors 

Previously,  pyrrolopyrimidine  core-containing  compounds  were 
reported  to  inhibit  HER  tyrosine  kinase  activity  [53,54].  Struc¬ 
tures  and  HER2  tyrosine  kinase  inhibitory  activities  of  some 
representative  compounds  are  shown  in  Figure  12  and  Table"^ 
[122].  Most  of  these  compounds  inhibited  HERl  and  -2  phos¬ 
phorylation  with  IGjq  values  <100  nM.  In  addition,  several  of 
these  compounds  also  inhibit  phosphorylation  of  VEGFR  fam¬ 
ily  members  with  IGjg  values  <  300  nM.  Hence,  they  are  consid¬ 
ered  as  dual  inhibitors  of  both  HER  and  the  VEGFR  family  of 
receptor  tyrosine  kinases.  Additional  analogs  of  pyrrolopyrimi- 
dines  as  inhibitors  of  HER,  VEGFR,  Bcr-Abl,  c-Abl  and  c-Raf-1 
kinases  were  also  reported  [123].  Three  representative  analogs  are 
shown  in  Figure  13.  Initially,  pyrrolopyrimidines  were  developed 
as  HERl  tyrosine  kinase  inhibitors.  Nevertheless,  now  they  are 
considered  as  dual  HERl  and  -2  inhibitors  and  also  show  con¬ 
siderable  inhibition  of  phosphorylation  of  VEGFR  and  other 
kinases  [54].  Q 

PK-II66  (98)  is  an  example  of  the  pyrrolopyrimidine  HER 
tyrosine  kinase  inhibitors  under  advanced  stages  of  clinical 
studies  (Figure  12)  [55].  PKII66  inhibits  kinase  activity  of 
HERl  and  -2  with  ICjq  values  of  1  and  1 1  nM,  respectively. 
PKI-166  inhibits  EGF  activated  HERl  phosphorylation  and 
activated  HERl -mediated  signalling  with  an  IC^q  value  of 
51  nM  [54].  PKI-166  also  inhibited  HER2  autophosphorylation 


in  two  hum^ 
with  IG, 


reast  cancer  cell  lines  (SK-BR-3  and  BT-474) 
of  0.1  -  1  pM.  In  cell-based  assays,  PKI-166 
inhibited  ^'^th  of  HERl -overexpressing  A431  and 
NCI-H596‘Cell  lines  and  HER2-overexpresing  SK-BR-3  and 
BT-474  cell  lines.  PKI-166  also  blocks  the  growth  of  prostate 
cancer  in  mice  in  an  androgen-dependent  manner  [56].  In  nude 
mice  implanted  with  human  pancreatic  carcinoma  cells, 
PKI-166  treatment  reduced  tumor  volume  by  45%,  and  in 
combination  with  gemcitabine  reduced  tumor  volume  by  85%, 
whereas  gemcitabine  alone  reduced  tumor  volume  by  only  59% 
[44].  In  a  recently  completed  Phase  I  study,  PKI-166  was  orally 
given  to  patients  with  advanced  solid  tumors.  PKI-166  showed 
acceptable  toxicity,  safety  and  pharmacokinetic  profiles. 
Dose-limiting  transaminase  elevations,  diarrhea,  skin  rash,  nau¬ 
sea  and  vomiting  were  the  principal  toxicides  of  PKI-166.  Fur- 
thermoire,  the  pharmacokinetic  profile  revealed  that  PKI-166  is 
orally  bioavailable,  quickly  absorbed  and  showed  a  linear 
dose-response  relationship  without  drug  accumulation  after 
multiple  doses.  Although  there  are  no  tumor  responses 
observed  with  PKI-166  treatment,  stable  disease  was  observed 
in  patients.  Additional  details  on  preclinical  and  clinical 
performance  of  PKI-166  are  awaited. 

2.3  Small-molecules  targeting  HER2  expression 

HER2  receptor  overexpression  is  associated  with  a  variety  of 
human  cancers.  HER2  overexpression  is  also  connected  to  the 
aggressive  nature  of  the  disease  and  short  survival  times.  Clin¬ 
ically  approved  HER2-specific  anticancer  therapeutics  are  tar¬ 
geted  at  elevated  levels  of  the  receptor.  It  seems  that  enhanced 
therapeutic  benefits  could  be  achieved  through  downregula- 
tion  of  HER2  by  inhibiting  HER2  gene  expression  in  addi¬ 
tion  to  targeting  elevated  levels  of  HER2  receptor  that  are 
already  expressed.  It  has  been  demonstrated  that  HER2  over¬ 
expression  is  accomplished  through  gene  duplication  and 
increased  transcription.  Ets  factor  ESX  (epithelial-restricted 
with  serine  box)  is  a  transcription  factor  that  activates  the 
HER2  gene  in  breast  cancers  [57].  ESX  is  overexpressed  in 
HER2-overexpressing  breast  cancers.  It  has  been  demon¬ 
strated  that  the  interaction  between  ESX  and  another  nuclear 
protein,  Sur-2/DRIP130  (a  Ras-linked  subunit  of  the  human 
mediator  complex),  is  required  for  overexpression  of  HER2  in 
malignant  breast  cancers.  Efficient  inhibition  of  this  interac¬ 
tion  reduces  HER2  expression.  Recently,  a  number  of  ada- 
manolol  derivatives  were  claimed  as  inhibitors  of  Sur-2 
[124,125].  Several  representative  analogs  are  shown  in  Figure  14. 
These  peptidomimetic  small  molecules  were  claimed  to 
inhibit  expression  of  the  HER2  receptor  through  disrupting 
ESX-Sur-2  interactions.  Several  of  these  compounds  inhib¬ 
ited  ESX-Sur-2  interactions  at  ICjq  values  in  the  range  of  3  - 
100  pM.  Limited  structure-  activity  relationships  on  some  of 
the  adamanolol  analogs  are  also  described. 

2.4  Small-molecule  inhibitors  of  HER2  shedding 

HER2  undergoes  a  slow  proteolytic  shedding  of  its  extracellu¬ 
lar  domain  in  cultured  cells.  The  soluble  HER2  extracellular 
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Figure  12.  Representative  pyrrolopyrimidine-based  HER2  kinase  inhibitors. 


Figure  13.  Representaive  analogs  of  pyrrolopyrimidine-based  FIER2  kinase  inhibitors. 


domain  was  found  in  serum  of  patients  with  advanced  breast 
cancer,  where  it  correlates  with  a  decreased  response  to  ther¬ 
apy.  In  HER2-expressing  clinically  poor  responsive  cancers, 
the  receptor  is  processed  by  an  unknown  matrix  metallopro- 
tease  to  yield  a  truncated  membrane  associated  receptor  (also 
known  as  p95HER2)  and  a  soluble  extracellular  domain 
[58-60].  Interestingly,  like  other  HER  family  receptors,  particu¬ 
larly  HERl,  loss  of  the  ligand  binding  extracellular  domain 
renders  the  HER2  intracellular  membrane  associated  domain 
a  constitutively  active  tyrosine  kinase.  Furthermore,  the 


increased  level  of  p95HER2  in  HER2-expressing  breast 
cancers  is  correlated  with  a  reduced  responsiveness  to  treat¬ 
ment  and  a  decreased  overall  clinical  outcome  [6i].  On  the 
basis  of  these  observations  it  has  been  hypothesized  that  the 
processing  of  the  HER2  extracelluar  domain  generates  a 
constitutively  active  receptor  kinase  that  can  continuously 
deliver  growth  and  survival  signals  to  cancer  cells.  A 
method  to  disrupt  shedding  of  HER2  extracellular  domain 
using  matrix  metalloprotease  antagonists  has  been  claimed 
[126].  In  a  proof  of  concept  study  using  GM-6001  (107),  a 
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Figure  14.  Reprsentative  peptidomimetic  small-molecule  Inhibitors  of  HER2  expression. 


Figure  15.  Structre  of  GM6001,  a  broad  spectrum  matrix 
metalloprotease  antagonist  that  inhibits  FIER2  shedding. 


matrix  metalloprotease  antagonist,  particularly  a  MMP-15 
antagonist,  reduced  shed  HER2  extracellular  domain  serum  lev¬ 
els  (Figure  15).  A  diagnostic  method  to  assess  elevated  level  of 
MMP-15  in  patients  with  HER2-expressing  cancer  is  also 
claimed  [126].  The  elevated  level  of  MMP-15  in  HER2-express- 
ing  cancer  indicates  that  the  patient  has  an  elevated  p95HER2  or 
shed  HER2  extracellular  serum  level  or  will  have  a  poor  clinical 
outcome.  The  small-molecule  compounds  that  disrupt  shedding 
of  the  HER2  extracellular  domain  and  generation  of  a  constitu- 
tively  active  intracellular  receptor  kinase  represent  a  novel  class  of 
anticancer  therapeutics  targeting  HER2-expressing  cancers. 


Table  7.  Inhibition  of  tyrosine  kinase  activity  oMER2 
and  HER1  by  the  pyrrolopyrimidine  analogs  90ac  97. 


Compound  IC50  (pM)  O 


HER2 

HER1  ^ 

90 

0.008 

0.0031*^ 

91 

0.0072 

0.003Q 

92 

0.0067 

0.0031 

93 

0.005 

0.007 

94 

0.0094 

O 

0.0024 

95 

0.005 

0.0043 

96 

0.005 

0.0047 

97 

0.0085 

0.0063 

3.  Expert  opinion 

The  HER  family  of  receptors  are  validated  anticancer  drug 
targets.  Profound  clinical  performance  of  well-designed 
first-generation  small-molecule  and  humanized  mAh  thera¬ 
peutics  signifies  the  rationale  behind  targeting  the  HER  fam¬ 
ily  of  receptors.  Several  investigational  drugs  targeting 
members  of  the  HER  family  are  in  various  stages  of  clinical 
development.  This  indicates  the  intensiveness  of  drug  discov¬ 
ery  efforts  targeting  HER  family  receptors  that  mediate  com¬ 
plex  and  vital  signalling  pathways.  The  long-term 
consequences  of  disruption  of  signalling  pathways  mediated 
by  receptor  tyrosine  kinases  using  small-molecule  or  human¬ 
ized  mAbs  have  yet  to  be  completely  understood.  For  exam¬ 
ple,  a  recent  report  on  cardiotoxic  effects  of  imatinib  mesylate. 
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the  first  Bcr-Abl  kinase  inhibitor  to  be  approved  for  use  in  the 
clinic,  highlight  the  surprising  outcomes  that  may  arise  from 
disruption  of  kinase  mediated  signalling  mechanisms  [62], 
Growing  understanding  of  the  structural  details  of  HER  fam¬ 
ily  receptor  kinase  domains  will  allow  the  design  of  highly 
selective  inhibitors  targeting  either  active  or  inactive  confor¬ 
mations  of  individual  receptors.  Finally,  considering  the 
prominent  role  of  the  HER2  signalling  pathway  in  cancer. 


inhibitors  of  expression  and  shedding  of  HER2  will  also  pro¬ 
vide  an  opportunity  to  study  the  consequences  of  selective 
inhibition  of  HER2. 
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Abstract 

Recent  studies  indicate  that  p53-MDM2  interactions  can  be  successfully  inhibited  by  small 
molecules  or  peptides.  In  this  work,  we  performed  pharmacophore  model  studies  including  the 
substructure  concept  and  the  explicit  shape  constrains  to  the  active  conformation  of  Nutlin  or  the 
p53  fragment  that  binds  to  the  MDM2.  Database  screening  using  these  models  resulted  in  a 
series  of  structurally  novel  compounds.  Eight  eight  representative  hits  were  selected  for  initial 
cytotoxicity  assay  against  HCT116  p53'^^'*‘  and  p53'^'  cell  lines.  Eight  of  the  compounds  showed 
desirable  activities  in  these  cells  as  well  two  additional  cell  lines. 


Keywords 
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Introduction 


Mouse  double  minute  2  (MDM2)  is  an  oncogene  that  was  originally  isolated  from  a 
spontaneously  transformed  mouse  fibroblast  cell  line/  Later  this  gene  was  rediscovered  as  a  p53 
binding  protein  in  several  rat  fibroblast  cell  lines.  The  human  homologue  has  been  shown  to  be 
over-expressed  in  40-60%  of  human  osteogenic  sarcomas  and  about  30%  of  soft  tissue 
sarcomas. The  human  MDM2  gene  encodes  a  491  amino  acid  residues  that  includes  a  p53- 
binding  domain,  an  acidic  region,  a  zinc-finger,  and  a  ring-finger  domain  (see  Figure  1)."^ 

p53  is  a  393  amino  acid  residue  transcription  factor  that  up-regulates  the  expression  of  proteins 
that  mediate  distinct  cellular  responses  such  as  apoptosis  and  cell  growth  arrest.^’^  Over  50%  of 
human  tumors  contain  mutations  in  p53  that  inactivate  its  transcription  functions.  Most  of  the 
p53  mutations  occur  within  the  DNA-binding  domain  such  that  p53  is  unable  to  act  as  a 

o 

sequence-specific  DNA  binding  transcription  factor.  The  fact  that  mutational  inactivation  of  p53 
was  observed  as  one  of  the  most  common  genetic  events  indicates  that  p53  plays  a  critical  role  in 
suppressing  the  tumor  growth.  Structural  studies  have  shown  that  p53  is  made  up  of  at  least  four 
functional  domains  regulating  its  transcriptional  activities  (see  Figure.  1):  1)  an  N-terminal  trans¬ 
activation  domain,  which  is  prolin-rich  and  required  for  interaction  with  components  of  the 
transcriptional  machinery,  2)  a  central  conserved  DNA-binding  core  domain,  3)  a 
tetramerization  domain  that  assists  sequence-specific  DNA  binding,  and  4)  a  C-terminial 
negative  regulatory  domain  whose  phosphorylation  primes  the  latent  sequence-specific  DNA- 
binding  function  of  p53  for  activation."^  The  N-terminal  transactivation  domain  contains  the 
binding  pocket  of  the  most  important  negative  cellular  regulators  of  the  p53,  i.e.  the  product  of 
the  MDM2  oncogene.^ 

The  MDM2  protein  regulates  the  p53  protein  activity  in  at  least  three  different  ways.^°  Briefly,  1) 

it  binds  to  the  p53  trans activation  domain,  blocking  the  transcriptional  activity,  2)  upon  binding 

to  p53,  it  induces  nuclear  export,  and  3)  it  stimulates  p53  degradation  by  catalyzing  the 
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ubiquitination  of  p53  through  the  uniquitin/proteasome  pathway,  as  reviewed  recently.  ’ 
Transcription  of  the  MDM2  gene  itself  is  regulated  by  p53  in  response  to  the  cellular  stress. 


indicating  p53  and  MDM2  form  an  autoregulartory  feedback  loop/°  Therefore,  as  one  of  the 
approaches  to  stabilize  the  p53  level,  MDM2  has  been  selected  as  a  very  attractive  therapeutic 
cancer  target.  Efforts  are  underway  to  design  peptides  or  small-molecule  compounds  to  disrupt 
the  p53-MDM2  interactions.  Biological  and  functional  studies  have  mapped  the  residues  of  19- 
102  as  the  domain  for  p53  binding,  while  MDM2  binding  domain  of  p53  was  narrowed  down  to 
15-residue  peptide."^  The  first  crystal  complex  of  MDM2  with  15-mer  wildtype  p53  was  resolved 
in  1996  and  clearly  revealed  that  MDM2  possesses  a  well-defined  deep  hydrophobic  binding 
pocket  efficiently  occupied  by  the  triad  motif,  i.e.  F19,  W23  and  L26,  of  p53.^^  In  addition  to  the 
H-bond  interaction  between  N  atom  of  W23  side  chain  of  p53  and  target  MDM2,  multiple  van 
der  Walls  contacts  were  observed  to  stabilize  the  complex. 

Most  of  the  data  on  the  inhibition  of  the  p53-MDM2  interaction  have  been  obtained  with 
peptides.  The  initial  work  suggested  that  p53  peptide  could  be  used  as  probes  to  investigate  the 
interaction  and  confirmed  the  critical  function  of  the  triad  motif  as  revealed  from  the  resolved 
p53-MDM2  complex.^’'"^’^^  Subsequently,  a  potent  12-mer  peptide  was  identified  by  screening 
phage  peptide  libraries. Several  shorter  peptides  showed  tighter  binding  as  observed  from  the 
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kinetic  and  thermodynamic  studies. 

Several  small  molecule  inhibitors  of  the  p53-MDM2  interactions  including  chalcones, 
chlorofusin,  etc,  have  also  been  discovered  recently. A  major  development  in  the  field  was 
the  discovery  of  Nutlin  series  of  compounds  and  the  determiniation  of  co-crystal  structure  at  2.3 
A  resolution.  The  Nutlin  binding  mode  as  revealed  in  the  crystal  structure  mimics  the 
interactions  of  the  p53  peptide  to  a  high  degree.  One  bromophenyl  moiety  sitting  deeply  in  the 
W23  pocket,  the  other  bromophenyl  group  overlapping  the  L26  space,  and  the  hydrophobic  ethyl 
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ether  side  chain  directed  toward  the  F19  pocket. 

Taken  together,  all  these  studies  strongly  prove  that  protein-protein  interactions  can  be 
successfully  manipulated  with  small  molecules.  Therefore,  p53  pathway  can  be  effectively 
targeted  to  develop  novel  anti-cancer  drugs. 


In  this  work,  we  carried  out  extensive  pharmacophore  model  studies  using  the  active 
conformations  of  imidazoline  analogue  (Nutlin)  and  the  side-chain  of  P53  fragment  which  were 


resolved  in  the  crystallographic  structures  of  MDM2  complex.  The  derived  models  were  applied 
to  screen  our  in-house  small  molecule  database  and  the  representative  compounds  were  tested  in 
HCT116p53^'^  and  HCTl  16  pSS'^' cell-lines. 

Computational  Methods 

Substructure  Search.  Initially,  we  applied  a  substructure  concept  to  find  a  spiral  compounds 
with  extended  arms,  including  Nutlin-like  molecules.  The  Nutlin  compound  was  built  by 
Catalyst  and  converted  to  a  hypothesis.  Three  additional  structural  motifs  mimicking  the  Nutlin 
central  bridge  were  designed,  and  these  hypothesis  were  used  as  queries  to  search  our  database  of 
small  molecule  compounds. 

Functional  Feature  Pharmacophore  Model.  Functional  feature  pharmacophore  model  studies 
were  extensively  performed  on  both  Nutlin  and  the  p53  side  chain  of  the  conserved  triad  motif, 
i.e.  F19,  W23,  L26  by  Catalyst  software  package.  The  conformation  of  Nutlin  or  p53  was  taken 
from  the  observation  as  resolved  in  the  x-ray  structures,  as  referred  to  as  active  conformation. 

In  light  of  Nutlin  pharmacophore  model  development,  we  first  mapped  H-bond  donor,  H-bond 
acceptor,  hydrophobic  and  aromatic  ring  features  to  the  active  conformation.  We  then  assigned 
geometrical  constraints  to  each  feature,  i.e.,  coordinate  and  size  of  the  feature.  Finally,  all  the 
selected  features  were  merged  into  one  pharmacophore  model.  Following  the  similar  procedure, 
pharmacophore  model  of  p53  triad  motif,  F19,  W23,  L26,  was  derived. 

Shape-merged  single  compound  pharmacophore  query.  The  Catalyst  was  used  to  generate 
the  explicit  shape  constraint  of  Nutlin  or  p53  side  chain  ligand.  Fist,  the  active  conformation  of 
the  template  molecule  (Nutlin  or  p53  fragment  of  amino  acid  residues  from  19-26)  was  mapped 
to  the  corresponding  feature  model,  respectively.  Then  the  explicit  shape  constraint  was 
generated.  Finally,  the  shape  was  merged  with  the  feature  model  as  a  single  query  in  Catalyst, 
that  we  referred  to  as  the  shape-merged  pharmacophore  model.  The  similarity  tolerance  of  shape 
is  an  important  parameter  to  adjust  the  percentage  of  hits  yielded  from  database  screening.  The 
values  of  0.5-1  was  used  as  default  value  in  the  shape-merged  p53  side  chain  model,  but  for  the 
shape-merged  Nutlin  feature  model,  we  relaxed  the  shape  similarity  to  values  of  0.45-1. 


Docking  studies.  The  docking  studies  were  carried  out  by  the  GOLD  (version  1.2  Genetic 
Optimization  for  Ligand  Docking)  software  package  running  on  our  multi-processor  linux  PC 
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and  a  24-processor  Silicon  Graphics  Onyx  workstation  as  described.  We  used  the  crystal 
structure  of  Nutlin  bound  MDM2  complex  (pdb  code  IRVl)  as  the  template  target,  and  the 
ligand,  Nutlin,  was  subsequently  removed  to  keep  the  binding  pocket  available.  Hydrogen  atoms 
were  added  to  both  protein  and  ligand  according  to  the  protonation  state  at  pH  7.0  during  the 
docking  studies. 
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The  docking  area  was  defined  as  a  sphere  with  the  radius  of  15  A  at  the  central  coordinate 
position  (  47.47,  12.31,  35.24)  of  the  Nutlin  molecule  as  chrystallographically  observed  in  the 
MDM2  binding  pocket.  The  ligand  was  placed  within  the  active  site  using  a  least- squares  fitting 
procedure.  All  docking  runs  were  performed  using  standard  default  settings  with  a  population 
size  of  100,  a  maximum  number  of  100,000  operations,  a  mutation  and  crossover  rate  of  95.  The 
scoring  function  was  contributed  basically  from  H-bond  interaction,  van  der  Waals  interactions 
within  the  complex,  and  the  ligand  internal  energy  which  was  summarized  by  ligand  steric  and 
torsional  energies.  Explicit  electrostatic  interaction  was  ignored  but  was  modeled  into  the  H- 
bond  interactions. 

Database  search.  Substructure  queries,  and  the  shape-merged  pharmacophore  models  derived 
from  either  the  Nutlin  or  the  side  chain  of  p53  fragment  in  their  active  conformation, 
respectively,  were  applied  to  screen  our  in-house  database  by  Catalyst. 

Cell.  Cancer  cells  (HCTl  16  p53^^'^,  HCTl  16  p53'^')  and  human  breast  cancer  cell  MDA-MB-435 
(p53  mutant)  were  obtained  from  Dr.  Bert  Vogelstein  (Johns  Hopkins).  The  human  ovarian 
carcinoma  cell  line  (HEY  p53  wildtype)  naturally  resistant  to  cisplatin  was  kindly  provided  by 
Dr.  Dubeau  (University  of  Southern  California  Norris  Cancer  Center).  Cells  were  maintained  as 
monolayer  cultures  in  RPMl  1640  supplemented  with  10%  fetal  bovine  serum  (Gemini- 
Bioproducts,  Woodland,  CA)  and  2  mmol/E  L-glutamine  at  37°C  in  a  humidified  atmosphere  of 
5%  CO2.  To  remove  the  adherent  cells  from  the  flask  for  passaging  and  counting,  cells  were 
washed  with  PBS  without  calcium  or  magnesium,  incubated  with  a  small  volume  of  0.25% 
trypsin-EDTA  solution  (Sigma-Aldrich,  St.  Eouis,  MO)  for  5  to  10  minutes,  and  washed  with 


culture  medium  and  centrifuged.  All  experiments  were  carried  out  using  cells  in  exponential  cell 
growth. 

Compounds.  The  selected  compounds  were  prepared  in  100%  DMSO  at  a  concentration  of 
lOmM  and  stored  at  -20°C.  Further  dilutions  of  each  drug  were  freshly  made  in  phosphate  saline 
buffer  to  final  concentration  before  use. 

Cytotoxicity  Assay.  Cytotoxicity  was  assessed  by  a  3-(4,5-dimethylthiazol-2-yl)-2,5- 
diphenyltetrazolium  bromide  assay  as  previously  described.  Briefly,  cells  were  seeded  in  96- 
well  micro-titer  plates  (4,000  cells/well)  and  allowed  to  attach.  Cells  were  subsequently  treated 
with  a  continuous  exposure  to  the  corresponding  drug  for  72  hours.  A  3-(4,5-dimethylthiazol-2- 
yl)-2,5-diphenyltetrazolium  bromide  solution  (at  a  final  concentration  of  0.5  mg/ml)  was  added 
to  each  well  and  cells  were  inculated  for  4  hours  at  37  °C.  After  removal  of  the  medium,  DMSO 
was  added  and  the  absorbance  was  read  at  570  nm.  All  assays  were  performed  at  least  in 
triplicate.  IC50  (dose  of  each  drug  that  inhibits  50%  of  cell  growth)  was  then  determined  from  a 
plot  of  log  (drug  concentration)  versus  percentage  of  cell  kill. 

Results  and  Discussions 

Substructure  Query.  Substructure  query  approach  of  active  compounds  is  an  easy  first  step  to 
rapidly  identify  hits.  The  recent  crystal  structure  (IRVl)  revealed  a  deep  substrate  binding 
pocket,  suggesting  that  a  variety  of  compounds  with  a  central  cyclic  structure  and  finger-like  side 
chain  substituents  could  occupy  the  cavity.  Therefore,  our  first  5-memember  ring  hypothesis,  sql 
(Table  1)  contained  the  central  core  of  Nutlin,  while  other  three  were  conceptually  bridge-like 
fragments  functioning  as  searching  queries  (sq2-sq4.  Table  1).  The  number  of  hits  from  two 
databases,  Chembrdige  and  Chemdiv,  against  each  substructure  motif  are  listed  in  Table  1. 

Pharmacophore  Model  Studies 

Nutlin  -based  pharmacophore  model.  Referring  to  the  interactions  between  Nutlin  and  the 
target  protein,  MDM2,  determined  in  the  X-ray  structure(lRVl)  ,  only  hydrophobic  features 
(P1-P6)  were  considered  in  the  pharmacophore  model  (  Figure  2a).  The  features  PI,  P2,  and  P4 
were  mapped  by  the  benzene  ring  and  P3  by  aliphatic  methane  group  which  was  observed  to 


interact  with  MDM2  by  its  Van  der  Waals  contacts.  Additionally,  P5  and  P6  were  recognized  by 
the  bromine  substitutes.  Shape-merged  pharmacophore  model  was  recently  studied  on  several 
potent  HIV-1  integrase  inhibitors  and  successfully  identify  novel  compounds  to  block  the 
integrase  activities.  We  applied  the  similar  concept  to  generate  the  explicit  shape  of  Nutline  in  its 
active  orientation,  as  shown  in  Figure  2b.  Apparently,  both  feature  model  (Figure  la)  and  shape 
constraint  itself  (Figure  2b)  can  be  used  independently  as  searching  query  to  screen  small 
molecule  database.  In  this  work  however,  we  combine  the  shape  constraint  onto  the  feature 
model  to  generated  shape-merged  pharmacophore  model,  labeled  as  SMI  (Figure  2b),  in  order  to 
comprehensively  extract  the  chemical  characteristics  of  the  potent  imidazole  antagonist. 
Consequently,  the  shaped-merged  model,  SMI,  was  applied  to  screen  our  in-house  database,  705 
compounds  were  identified.  We  selected  the  compounds  with  wide  structure  variant  for  cell  line 
analysis.  The  structures  and  the  predicted  docking  score  are  listed  in  Table  2  and  some  of 
cytotoxic  measurement  of  the  selected  compounds  are  listed  in  Table  4. 

p53  side  chain  -based  pharmacophore  model 

Since  the  p53-MDM2  fragment  was  resolved  by  x-ray  technique,  a  large  number  of  peptide 
library  was  designed  to  compete  the  p53  binding.  The  first  published  putative  p53-based 
pharmacophore  model  addressed  the  hydrophobic  features  recognized  by  the  side  chain  of  F19, 
W23  and  L26,  and  H-bond  donor  mapped  by  the  W23  indole  side  chain.  Its  application  to 
screen  NCI  database  resulted  series  of  non-peptide  sulfonamide  inhibitors.  In  this  work, 
however,  we  extended  the  pharmacophore  studies  by  taking  into  account  explicit  shape  of  the 
active  conformation  of  the  p53  fragment,  determined  by  the  x-ray  structure,  therefore,  two  shape- 
merged  pharmacophore  models  were  systematically  generated,  as  shown  in  Figure  3.  Figure  3a 
shows  the  model  with  the  consideration  of  only  three  hydrophobic  features,  but  enhanced  by 
shape  constraint,  referred  to  as  SM2.  Figure  3b  shows  the  shape-merged  model,  referred  to  as 
SM3,  that  keeps  additional  hydrogen-bond  feature.  As  consequence,  SM2  and  SM3  were  treated 
as  independent  queries  to  virtually  screen  our  in-house  database.  The  structure  and  GOLD 
predicted  docking  scores  of  selected  molecules  identified  by  the  pharmacophore  models  are 
listed  in  Table  3,  and  the  cytotoxicity  measurements  are  in  Table  4. 


Docking  Studies 

Docking  simulations  were  performed  on  several  potent  inhibitors,  to  predict  the  binding 
mechanism  (Figure  4) 


Conclusions 
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Table  1  Substructure  motif  used  as  database  searching  query  and  the  number  of  hits 


Table  2  Structures  and  GOLD  docking  scores  of  selected  compounds  identified  from 


substructure  queries. 
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Table  3  Structures  and  GOLD  docking  scores  of  selected  compounds  identified  from 


Nutlin  and  p53  side  chain  pharmacophore  models 
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Table  4  Cytotoxicity  assay  of  selected  compounds  in  a  panel  of  cancer  cell  lines 
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Figure  1 .  Schematic  diagram  showing  the  functional  domains  of  p53  and  MDM2. 


Figure  2.  Pharmacophore  models  generated  from  active  conformation  of  Nutlin.  (a) 
Mapping  Nutlin  onto  the  feature  pharmacophore  model,  (b)  Mapping  Nutlin  to  the 
shape-only  query,  (c)  Mapping  Nutlin  onto  shape-merged  pharmacophore  model. 
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Figure  3.  Pharmacophore  models  derived  from  the  side  chains  of  p53 
fragment  interacting  with  MDM2  p53  binding  domain,  (a)  p53  fragment 
(amino  acide  residue  19-26)  mapping  onto  the  3-hydrophobic  feature 
pharmacophore  model,  (b)  p53  fragment  (amino  acide  residue  19-26) 
mapping  onto  the  hetero-feature  pharmacophore  model. 
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Figure  4.  Potent  molecules  binding  mode  analysis  predicted  from  docking  study.  MDM2  binding  surface 
is  displayed  in  yellow  (dot  or  solid)  surface,  and  ligand  is  rendered  as  stick  or  ball-stick,  (a)  and  (b) 
show  the  MP15  binding  into  the  MDM2  pocket,  (a  is  for  top  view  and  b  for  side  view),  (c)  and  (d)  show 
the  MP17  binding  into  the  MDM2  pocket,  (e)  and  (f)  show  the  MP22  binding  into  the  MDM2  pocket,  (g) 
and  (h)  show  the  MP31  binding  into  the  MDM2  pocket. 
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Introduction 


The  human  epidermal  growth  factor  receptor- 2  (HER2)  proto-oncogene  encodes  a  185- 
kDa  glycoprotein  that  plays  a  key  role  in  breast  cancer.  The  HER  tyrosine  kinase  receptor 
family  comprises  four  homologous  epidermal  growth  factor  (EGE)  receptors: 
EGER/ErbBl/HERl,  HER2/Neu/ErbB2,  HER3/ErbB3,  and  HER4/ErbB4.  The  HER 
family  receptors  are  among  the  most  studied  cell  signaling  families  in  cancer  biology.' 
Each  protein  has  an  extracellular  ligand-binding  domain,  a  transmembrane  region,  and  an 
intracellular  cytoplasmic  domain.  The  extracellular  domain  is  involved  in  recognizing 
and  binding  ligands  that  are  able  to  activate  receptor.  The  cytoplasmic  domain  contains  a 
tyrosine  kinase  region  and  a  carboxy  terminal  tail  that  harbors  tyrosine 
autophosphorylation  sites.  The  HER  family  receptors  have  similar  structures  and  a  high 
degree  of  homology  in  the  tyrosine  kinase  region,  although  they  possess  diverse 
characteristics  that  determine  their  signaling  specificity.  The  tyrosine  kinase  domains  of 
HER2  and  HER4  show  approximately  80%  homology  to  that  of  HERl,  whereas  HER3 
lacks  intrinsic  tyrosine  kinase  activity.  Interestingly,  the  extracellular  domains  are  less 
conserved  among  the  four  receptors,  which  is  indicative  of  different  specificity  in  ligand 
binding. 

The  HER  receptors  exist  as  inactive  monomers.  A  variety  of  ligands  can  bind  to 
extracellular  domains  of  HERl,  HER3  and  HER4  and  activate  these  receptors.  These 
ligands  have  different  specificities  for  each  receptor,  resulting  in  different  cellular 
signaling  effects."'  There  is  no  known  ligand  for  HER2,  rendering  it  as  an  orphan  receptor. 
Upon  binding  to  ligands,  the  extracellular  domains  of  HERl,  HER3,  and  HER4  receptors 


undergo  conformational  changes  that  allow  activation  of  the  receptors  by  either  a 
homodimerization  or  a  heterdimerization  process,  which  is  essential  to  stimulate  the 
tyrosine  kinase  activity  of  the  receptors  for  subsequent  generation  of  an  intracellular 
signal.  The  pairing  of  two  receptors  of  the  same  type  results  in  ‘homodimers’,  whereas 
pairing  of  two  different  type  of  receptors  produces  ‘heterodimers’.  Interestingly, 
constitutive  conformation  of  the  extracellular  domain  of  HER2,  an  orphan  receptor, 
resembles  the  conformation  of  ligand  activated  extracellular  domains  of  the  other  HER 
receptors.  This  autoactivated  conformation  enables  HER2  to  perform  as  a  preferred 
dimerization  partner  for  all  the  other  HER  receptors.^’  ^ 

In  a  variety  of  human  cancer  cells,  aberrant  signaling  involving  HER  receptors 
stimulate  pathways  that  activate  many  of  the  properties  associated  with  cancer,  including 
proliferation,  migration,  metastasis,  angiogenesis,  and  resistance  to  apoptosis.  Due  to  the 
high  frequency  of  abnormalities  in  receptor  signaling  in  human  cancers,  the  HER  family 
is  an  attractive  target  for  therapeutic  development.  The  frequently  observed  causes  of 
signaling  abnormalities  in  cancers  involving  HER  receptors  are  overexpression  of 
receptors,  overproduction  of  growth  factor  ligands,  and  ligand- independent  receptor 
activation.  The  HER  receptors  are  overexpressed  or  deregulated  in  a  wide  variety  of 
cancers  including  breast,  colorectal,  ovarian,  prostate,  and  non-small  cell  lung  cancers. 
Importantly,  overexpression  of  HER  receptors  is  associated  with  poor  disease  prognosis 
and  reduced  survival.  Mounting  pre-clinical  and  clinical  evidence  supports  the  rationale 
behind  the  HER  family  targeted  anticancer  therapeutic  approaches.  Although  the 
complexities  of  the  signaling  pathways  involving  the  HER  family  of  receptors  are  not 


fully  understood,  several  possible  points  within  the  pathways  have  been  identified  to  be 
therapeutically  advantageous  if  interrupted/  Drugs  targeting  HER  family  fall  into  three 
main  categories  depending  on  the  receptor  region  targeted:  extracellular,  intracellular  and 
nuclear.  The  most  promising  and  advanced  therapeutic  strategies  targeting  HER  receptors 
are  monoclonal  antibodies  and  small-molecule  tyrosine  kinase  inhibitors.  Intensive 
research  over  the  last  20  years  in  these  areas  has  yielded  a  number  of  promising 
anticancer  therapeutics.  Currently,  five  anticancer  therapeutics  targeting  the  HER 
receptor  family  have  been  approved.  Two  of  them  are  humanized  antibodies,  trastuzumab 
(Herceptin)  and  cetuximab  (Erbitux),  targeted  to  the  extracellular  domains  of  HER2  and 
HERl,  respectively.  Remaining  three  of  them  are  small- molecule  tyrosine  kinase 
inhibitors,  Gefitinib  (1),  and  Erlotinib  (2),  targeted  to  the  HERl  tyrosine  kinase  domain 

o 

and  GW2016/Tykerb/Eapatanib  (3)  targeted  to  the  HER2  tyrosine  kinase  domain 
(Figurel).  Several  small-molecule  HER  tyrosine  kinase  inhibitors  are  being  investigated 
in  various  stages  of  clinical  trials. These  compounds  compete  with  ATP  for  binding  to 
the  ATP  site  on  the  receptor’s  intracellular  tyrosine  kinase  domain.  By  binding  to  the 
ATP  site  and  inhibiting  catalytic  activity,  they  block  autophosphorylation  of  key  tyrosine 
residues  on  the  receptor  kinase  domain.  This  disrupts  activation  of  the  receptor  mediated 
downstream  signaling  cascades.  Considering  the  enormous  research  activity  in  the 
discovery  of  anticancer  therapeutics  targeting  the  HER  receptor  family  and  small 
molecule  inhibitors  that  target  the  intracellular  HER2  appear  to  be  a  very  promising 
approach  towards  treating  HER2  mediated  cancers.  These  molecules  act  by  binding  either 
reversibly  or  irreversibly  to  the  C-terminal  tyrosine  kinase  domain  of  HER2,  thereby 


inhibiting  autophosphorylation  of  the  receptor  and  therefore  activation.  Here,  the  study 
emphasis  on  HER2  targeted  small-molecule  therapeutics. 

In  the  drug  discovery  process,  computational  methodologies  are  well-established 
useful  tools.  Other  technologies,  such  as  combinatorial  chemistry  and  high  throughput 
screening  (HTS)  can  be  used  to  synthesize  and  test  thousands  of  compounds.  A  variety  of 
computational  techniques  are  currently  available  to  rapidly  screen  the  compound 
libraries.  Virtual  screening  uses  computational  models  to  predict  biological  activity  of 
compounds  from  existing  databases  or  virtual  libraries.  This  results  in  identification  of 
molecules  against  a  drug  target  of  interest.  There  are  two  basic  approaches  in  virtual 
screening  depending  on  the  availability  of  3D  structure  of  the  drug  target.  If  a  3D 
structure  of  the  drug  target  is  available,  either  from  an  experimental  (X-ray,  NMR)  or 
theoretical  (homology  modeling)  technique,  structure-guided  drug  design  procedures 
such  as  docking  can  be  applied.  Otherwise,  ligand-based  drug  design  methods  for 
example  QSAR  and  pharmacophore  models  can  be  used.  The  two  screening  approaches 
could  be  used  separately  or  combined."  In  this  context,  with  an  aim  to  discover  small- 
molecule  HER2  inhibitors,  we  employed  a  consensus  model  based  on  pharmacophore 
and  docking  simulations.  In  this  study  we  also  present  novel  HER2  inhibitors  with 
activity  in  SKBr3  cells  overexpressing  HER2. 


Materials  and  Methods 


Database  Mining  and  Ligand  Preparation.  Currently,  there  is  only  one  FDA  approved 
tyrosine  kinase  inhibitor  (Tykerb/GW2016)  (3)  selectively  targeting  HER2  available  for 
cancer  treatment.  Apart  from  this,  there  are  two  clinically  approved  HERl  specific 
tyrosine  kinase  inhibitors  gefitinib  (1),  and  erlotinib  (2)  also  that  inhibit  HER2 
phosphorylation  and  downstream  signaling  cascades  in  HER2/HER3  overexpressing 
cancer  cells.  It  has  been  observed  that  most  of  the  HER  tyrosine  kinase  inhibitors  in 
clinical  development  also  show  cross  receptor  activity.  Tykerb  (GW-2016)  (3)  is  a 
reversible  dual  HER2  and  HERl  tyrosine  kinase  inhibitor.  The  quinazoline  derivative 
represents  a  first  generation  HER2  targeted  small-molecule  anticancer  agent. 

Our  in  house  HER  database  consists  of  4500  HER2  antagonists  collected  from 
250  references.  1200  molecules  have  IC50  values  for  human  HER2.  In  order  to  develop 
a  virtual  screening  model  targeting  HER2  around  474  molecules  were  selected  based  on 
experimentally  determined  IC50  values,  structural  diversity  and  similar  assay 
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conditions. 

All  ligand  structures  were  built  in  Cerius2  and  minimized  using  Open  Eorce  Eield 
(OEE)  methods  with  the  steepest  descent  algorithm."^^  A  gradient  convergence  value  of 
0.001  k.cal/mol  was  used.  Partial  atomic  charges  were  calculated  using  the  Gasteiger 
method.  The  resulting  geometries  with  partial  charges  were  used  for  docking  calculations 
and  the  mol  files  of  all  molecules  were  exported  and  minimized  using  modified 
CHARMM  force  field  in  Catalyst  software  package. Conformational  analysis  was 


carried  out  with  Confirm  module.  Poling  algorithm  of  Confirm  module  was  used  to 
generate  250  conformations  for  each  molecule  with  in  20  k.cal/mol."^"^ 

Pharmacophore  Model  Generation  and  Validation.  Pharmacophore  modeling 
correlates  activities  with  the  spatial  arrangement  of  various  chemical  features.  A  set  of 
474  HER2  inhibitors  with  activity  data  (IC50)  spanning  over  5  orders  of  magnitude  (pICjo 
=  3-9)  were  selected.  The  data  set  were  divided  into  training  and  test  set.  The  training  set 
was  selected  by  considering  both  structural  diversity  and  wide  activity  range.  Highly 
active,  moderately  active  and  also  inactive  compounds  were  included  in  order  to  obtain 
critical  information  on  pharmacophore  requirements.  This  training  set  was  used  to 
generate  quantitative  pharmacophore  models  while  some  of  the  highly  active  molecules 
were  used  to  generate  qualitative  pharmacophore  models. 

Comman  feature  hypotheses  (qualitative  models)  were  produced  by  comparing  a 
set  of  conformational  models  and  a  number  of  3D  configurations  of  chemical  features 
shared  among  the  training  set  molecules.  This  results  in  a  qualitative  model  wherein 
important  chemical  features  can  easily  be  identified.  For  chemically  meaningful  patterns 
it  is  important  to  identify  such  chemical  features  before  proceeding  to  the  quantitative 
model  generation.  In  order  to  confirm  essential  features  among  HER2  inhibitors,  10 
common  feature  hypotheses  were  generated  using  most  active  molecules  (3-9)  (Figure  2). 
In  all  these  hypotheses  hydrogen  bond  acceptors  (HBA),  hydrogen  bond  donor  (HBD) 
and  hydrophobic  aromatic  (HRA)  groups  are  common  features.  However,  these  models 
cannot  be  directly  used  for  prediction  of  biological  activities  for  the  compounds  obtained 


from  database.  Therefore,  we  used  quantitative  model  that  can  predict  the  biological 
activities  of  compounds. 

In  the  quantitative  model,  the  22  compounds  training  set  (3,  7-27)  was  selected 
with  the  above  criteria  (Figure  1,  2  and  3).  A  database  of  452  known  HER2  inhibitor  was 
used  as  a  test  set  to  validate  quality  of  the  pharmacophore.  Compounds  were  classified 
as:  highly  active  (pICjo  >6.5)  moderately  active  (pICjo  =  6.5-5)  and  inactive  (pICjo  <5). 
While  generating  quantitative  model,  a  minimum  of  0  to  a  maximum  of  5  features 
involving  HBA,  HBD  and  HRA  features  were  selected  and  used  to  build  a  series  of 
hypotheses  using  default  uncertainty  value  of  3.  In  general,  pharmacophore  models 
should  be  statistically  significant,  predict  activity  of  the  molecules  accurately  and  identify 
active  compounds  from  a  database.  Therefore,  the  derived  pharmacophore  map  was 
validated  using  cost  analysis,  test  set  prediction,  enrichment  factor,  and  goodness  of  hit. 

Qualitative  pharmacophore  model  and  quantitative  pharmacophore  model  were 
developed  using  HipHop  and  HypoGen  modules  respectively,  implemented  in  the 
Catalyst  software  package. 

Modeling  of  Active  and  Inactive  States  of  HER2  Structure.  The  HER  receptors  show 
close  homology  in  their  sequence  and  domain  organization:  they  have  a  ~600-residue- 
long  extracellular  region  with  four  domains  responsible  for  ligand  binding,  a  single 
transmembrane  helix,  and  an  ~500-residue  long  intracellular  region  with  the  most 
conserved  tyrosine  kinase  domain,  followed  by  a  less  conserved  regulatory  tail  at  the  C- 


terminal  end  of  the  protein.  To  provide  a  moleeular  basis  for  understanding  the  role  of 
HER2  in  signalling  and  cancer,  we  tried  to  model  active  and  inactive  states  of  the 
receptor  based  on  co-crystal  structures  of  HERl  (1M17  PDB"^^and  IXKK  bound 

to  Erlotinib  (OSI-774)  (2)  and  Tykerb  (GW572016)  (3)  respectively.  The  overall  fold  of 
these  two  crystal  structures  is  similar,  but  there  are  significant  differences  in  the 
orientation  of  the  N  and  C-terminal  lobes,  the  C-terminal  tail,  and  helix.  These 
orientations  define  the  shape  of  the  ATP  binding  pocket  and  conformation  of  the 
activation  loop.  The  N  and  C-terminal  lobes  of  kinases  are  connected  by  a  flexible  hinge 
region.  The  relative  orientation  of  the  two  lobes  with  respect  to  each  other  influences  the 
shape  of  the  ATP  binding  cleft  and  is  dependent  on  the  activation  state  of  the  kinase  and 
the  presence  of  ligands.  The  ATP  binding  cleft  in  the  HERl-Tykerb  co-crystal  structure 
is  in  a  relatively  closed  conformation.  The  ATP  binding  cleft  in  the  HERl -Erlotinib  co¬ 
crystal  structure  is  in  a  more  open  conformation.  The  N-  terminal  lobe  in  HERl-Tykerb 
complexed  is  rotated  -12°  relative  to  its  position  in  HERl-Erlotinib.  The  magnitude  of 
the  difference  suggests  that  the  two  inhibitors  target  different  forms  of  the  enzyme,  active 
and  inactive.  In  HERl-Tykerb,  residues  971  to  980  of  the  C-  terminal  tail  form  a  short 
alpha-helix  that  packs  along  the  hinge  region  connecting  the  N-  and  C-terminal  lobes. 
This  helix  partially  blocks  the  front  of  the  ATP  binding  cleft.  A  second  helical  segment 
containing  residues  983  to  990  extends  along  the  NH2-  terminal  lobe  of  the  protein. 
HERl-Erlotinib  has  active  activation-loop  conformations  where  as  HERl-Tykerb  has 
inactive  activation-loop  conformations.'^^  Based  on  these  two  co-crystal  structures, 
1M17  and  IXKK,  we  tried  to  obtain  the  3D  model  of  HER2  in  active  and  inactive  states 
of  HER2.  The  amino  acid  sequence  of  HER2  kinase  domain  was  taken  from  Swissprot 


database  (P04626).  ClustalW  program  was  used  to  identify  the  homologous  regions 
between  the  two  proteins.  The  catalytic  domain  of  the  EGFR  and  HER  receptors  is  well 
conserved  with  80%  sequence  identity  (Figure  4).  Amino  acid  residues  Asp746  and 
Cys751  play  a  prominent  role  in  ATP  binding  pocket  and  confer  selectivity  on  EGFR  and 
HER2  ligands. 

Entire  modeling  was  done  using  Insight  II  suite  of  softwares  (Accelrys,  Inc). 
Modeler  algorithm  was  used  to  build  both  active  and  inactive-like  conformations  of 
HER2.  Modeler  is  an  automated  comparative  modeling  program  designed  to  find  the 
most  probable  structure  of  a  protein  sequence,  given  its  alignment  with  related  structures. 
The  model  was  obtained  by  the  optimal  satisfaction  of  spatial  restraints  derived  from  the 
alignment  and  was  expressed  as  a  probability  density  function  for  the  features  restrained. 
The  optimization  procedure  is  a  variable  target  function  method  that  applies  conjugate 
gradients  algorithm  to  position  all  non-hydrogen  atoms.  The  structure  has  been  further 
refined  using  molecular  dynamics  simulation  in  explicit  water  following  standard 
procedure.  Both  the  resultant  models  have  the  correct  stereochemistry  as  judged  by  the 
Ramachandran  plot  and  conserves  the  topological  and  active  site  features  of  the  active 
and  inactive-like  conformations  of  HERE  The  DEG  motif  of  kinases  is  a  part  of  the 
activation- loop  and  involved  in  coordination  of  ATP.  The  side  chain  conformations  of 
these  amino  acids  were  checked  and  found  to  be  in  similar  orientation  as  that  of  the 
template  structures.  The  homology  models  of  HER2  active  and  inactive  states  were  built 
to  explore  whether  inhibitors  favorably  bind  to  active  or  inactive  states  of  HER2. 


Moreover,  we  wanted  to  determine  which  state  of  HER2  (active  or  inactive)  is  suitable 
for  docking  simulations. 

Docking  Studies.  Docking  calculations  were  performed  using  GOLD,  Glide,  and  eHITS 
as  described  below  on  the  modeled  active  and  inactive  states  of  HER2.  Eor  each  ligand, 
the  best  20  scored  poses  were  saved. 

1.  GOLD."^^  A  collection  of  474  minimized  HER2  inhibitors  were  docked  on  to  HER2 

o 

active  and  inactive  state  with  the  GOLD  software  in  a  20A  radius  surrounding  the  active 
site  and  ranked  according  to  its  score.  Standard  set  parameters  of  GOLD  were  used 
through  out  simulations.  Lor  each  of  the  10  independent  genetic  algorithm  runs,  with  a 
selection  pressure  1.1,  100000  operations  were  performed  on  a  set  of  5  islands  with  a 
population  size  of  100  individuals.  Default  operator  weights  were  used  for  crossover, 

o 

mutation,  and  migration  of  95,  95  and  10  respectively.  Default  cut-offs  values  of  2.5  A 

o 

for  hydrogen  bonds  and  4.0  A  for  van  der  Waals  were  employed.  All  other  values  were 
set  to  the  default.  Top  20  poses  were  saved  for  each  ligand  and  best  score  values  were 
taken  to  correlate  with  experimental  data. 

2.  Glide."^^  Glide  XP  calculations  were  performed  with  Impact  version  3.5.  The  protein 
charged  groups  that  were  not  located  in  the  binding  pocket  nor  involved  in  salt  bridges 
were  neutralized.  The  center  of  the  grid  enclosing  a  box  was  defined  by  the  center  of  the 
bound  dummy  ligand  fragment  at  the  hinge  region  of  HER2  protein.  The  enclosing  box 

o 

and  bounding  box  dimensions  were  fixed  to  14  and  10  A,  respectively.  No  further 
modifications  were  applied  to  the  default  settings.  A  total  of  20  poses  were  saved  for  each 


ligand  and  calculated  Glide  score  for  conformations  of  each  molecule.  High  score  value 
for  each  ligand  was  selected  to  correlate  to  its  biological  activity  data. 

3.  eHITs.^°  Docking  was  performed  for  these  474  HER2  inhibitors  on  to  HER2  active 
and  inactive  states  by  using  the  eHiTS  docking  program.  eHITS  evaluates  all  the  possible 
protonation  states  for  the  receptor  and  ligands  automatically  for  every  receptor-ligand 
pair.  The  docking  method  systematically  covers  the  part  of  the  conformational  and 
positional  search  space  to  avoid  severe  steric  clashes.  The  top  20  conformations  were 
saved  and  eHITs  scores  were  calculated  for  each  of  the  20  saved  ligand  conformations. 
High  score  value  of  each  ligand  was  selected  to  correlate  to  its  biological  activity. 

Consensus  Model  Generation."  To  generate  a  consensus  model,  a  multiple  regression 
analysis  was  carried  out  using  pharmacophore  predicted  activity  (Hypo-5)  along  with  the 
three  different  scores  from  GOED,  Glide  and  eHITS  as  descriptors  (X))  and  experimental 
pICjo  as  a  dependant  variable  (Y).  Of  the  474  ligands,  158  were  used  in  the  training  set 
and  316  in  the  test  set.  The  combination  of  predicted  activity  from  pharmacophore, 
GOEDscore,  Glidescore  and  eHITSscore  gave  a  better  model  than  any  other  models  with 
independent  and/or  combination  of  scores  and  predicted  activity  from  pharmacophore 
(Hypo5  model). 

Database  Screening.  The  Hypo5  model  was  used  to  screen  a  database  of  350,000 
compounds.  The  search  retrieved  961  hits.  We  selected  531  compounds  with  a  predicted 
IC50  values  <1  pM,  for  docking  studies.  GOED,  Glide  and  eHITs  scores  were  calculated 
for  all  the  531  molecules  using  active  and  inactive  state  of  HER2.  All  these  values  were 


then  substituted  in  the  Eq.  1  and  Eq.  2  and  calculated  predicted  activity  value  for  each 
compound.  Einally,  57  active  hits  with  a  predicted  activity  values  of  <100  pM  were 
selected  for  in-vitro  assay  against  human  breast  cancer  cell  line,  SKBrS,  which  is  known 
to  overexpress  HER2.^''^^ 

Cell  Culture.  The  human  breast  cancer  cell  line  SKBrS  obtained  from  American  Type 
Culture  Collection  (Manassas,  VA)  was  used  for  the  cytotoxicity  assay.  Cells  were 
maintained  as  monolayer  cultures  in  RPMI  1640  supplemented  with  10%  fetal  bovine 
serum  (Gemini-Bioproducts,  Woodland,  CA)  and  2  mM  E-glutamine  at  37°C  in  a 
humidified  atmosphere  of  5%  CO2.  To  remove  the  adherent  cells  from  the  flask  for  pass 
aging  and  counting,  cells  were  washed  with  PBS  without  calcium  or  magnesium, 
incubated  with  a  small  volume  of  0.25%  trypsin-EDTA  solution  (Sigma-Aldrich,  St. 
Eouis,  MO)  for  5  to  10  minutes,  and  washed  with  culture  medium  and  centrifuged.  All 
experiments  were  done  using  cells  in  exponential  cell  growth.  A  10  mM  stock  solution  of 
each  compound  was  prepared  in  DMSO  and  stored  at  -80°C.  Eurther  dilutions  were 
freshly  made  in  media. 

Cytotoxicity  Assay.  Cytotoxicity  was  measured  using  a  3-(4,5-dimethylthiazol-2-yl)-2,5- 
diphenyltetrazolium  bromide  (MTT)  assay.  Briefly,  8,000  cells  were  seeded  in  96-well 
microtiter  plates  and  allowed  to  attach  over  night.  Cells  were  subsequently  treated  with  a 
continuous  exposure  to  the  corresponding  drugs  for  72  hours.  A  MTT  solution  (at  a  final 
concentration  of  0.3  mg/mL)  was  added  to  each  well  and  cells  were  incubated  for  4  hours 


at  37  °C.  After  removal  of  the  medium,  DMSO  was  added  and  the  absorbance  was  read  at 


570  nm.  The  IC50  values  were  then  determined  for  each  drug  from  a  plot  of  log  value  of 
drug  concentration  versus  percentage  inhibition  of  cells. 

Results  and  Discussion 

Pharmacophore  Models.  The  best  common  feature  pharmacophore  model  indicated  the 
importance  of  HBA,  HBD  and  HRA,  which  were  further  confirmed  in  the  quantitative 
models.  In  the  quantitative  model,  several  hypotheses  were  generated  using  the  most 
diverse  training  set  (3,  7  -  27)  along  with  experimental  activities  (Tablel).  Top  ten 
hypotheses  were  comprised  of  HBA,  HBD,  HRA  features.  The  statistical  values  of  the  10 
hypotheses  such  as  cost  values,  correlation  (r)  and  root  mean  square  deviations  (RMSD) 
are  given  in  Table  2.  In  addition  to  estimation  of  activity  of  training  set  molecules,  the 
pharmacophore  model  should  also  estimate  the  activity  of  the  test  set  compounds 
accurately.  Two  statistical  methods  were  employed  to  rank  the  ten  hypotheses.  In  the  first 
method,  all  the  10  hypotheses  were  evaluated  using  a  test  set  of  452  known  HER2 
inhibitors,  which  were  not  included  in  the  training  set.  These  molecules  are  covering  a 
wide  range  of  activities  spanning  from  pICso  =  3  to  9.  Predicted  activities  of  test  set 
were  calculated  using  all  ten  hypotheses  and  correlated  with  experimental  activity.  Out  of 
10  hypotheses,  Hypo5  showed  better  correlation  coefficient  than  other  nine.  A  second 
statistical  test,  including  calculation  of  false  positives,  false  negatives,  enrichment  and 
goodness  of  hit  were  calculated  to  determine  the  robustness  of  the  hypotheses.  Under  all 
validation  conditions,  Hypo5  showed  better  results  than  the  other  nine  hypotheses. 


On  the  basis  of  the  Hypo5  model,  predicted  IC50  values  for  22  training  set 
compounds  along  with  fitness  score  are  presented  in  Table  1.  In  the  training  set,  all  six 
highly  active  compounds  were  correctly  predicted  as  highly  active.  Of  the  seven 
moderately  active  compounds,  six  were  predicted  as  moderately  active,  while  one 
compound  was  predicted  as  highly  active  and  all  nine  inactive  compounds  were  predicted 
as  inactive.  From  Table  1,  it  was  also  noticed  that  compounds  were  not  only  correctly 
predicted  as  high,  moderate  and  low  activity  but  also  the  fit  values  confer  a  good  measure 
of  how  well  the  defined  features  of  Hypo5  fit  with  the  molecule.  In  the  training  set,  it  was 
observed  that  all  features  of  Hypo5  (HBA,  HBD  and  three  HRA)  were  mapped  on  to  the 
highly  active  compounds  but  in  moderately  active  compounds,  one  or  more  features  did 
not  map  well  or  were  partially  missing.  In  the  inactive  compounds  set,  one  or  more 
features  were  totally  missing  (Figure  5). 

Database  mining  was  performed  using  the  BEST  flexible  searching  technique. 
Hypo5  was  used  to  screen  a  database  of  known  HER2  inhibitors.  Based  on  the  data  and 
considering  a  compound  pICjo  >  6.5  as  a  highly  active  compound,  number  of  parameters 
such  as  hit  list  (Ht),  number  of  active  percent  of  yields  (%Y),  percent  ratio  of  actives  in 
the  hit  list  (%A),  enrichment  factor  (E),  false  negatives,  false  positives  and  goodness  of 
hit  score  (GH)  were  calculated  (Table  3).  Erom  the  table  3,  64  false  positives  and  55  false 
negatives  were  observed.  An  enrichment  factor  of  1.67  and  a  GH  score  of  0.536  indicate 
the  quality  of  the  model.  Scatter  plot  of  predicted  activity  obtained  from  Hypo5  against 
the  experimental  activity  (pICjo)  of  the  training  and  test  set  is  shown  in  Eigure  6.  Erom 
this  plot  it  was  noticed  that  the  Hypo5  model  has  a  greater  tendency  to  show  false 


positives.  To  avoid  false  positives  and  limit  false  negatives  and  also  to  understand  the 
binding  mechanism  of  inhibitors  to  HER2,  we  further  extended  this  study  to  structure 
based  ligand  design. 

Validation  of  HER2  Active  and  Inactive  States  and  Docking  Simulations.  We 

constructed  homology  models  for  HER2  active  and  inactive  state  conformations  to 
predict  the  correct  binding  mode.  Eirst,  the  3D  models  of  HER2  (active  and  inactive 
states)  were  validated  by  docking  the  Erlotinib  (the  bound  ligand  in  co-crystal  structure 
of  active  state  of  HERl  (1M17  PDB))  on  to  HER2  active  state  model  and  Tykerb  (the 
bound  ligand  in  co-crystal  structure  of  inactive  state  of  HERl  (IXKK  PDB))  on  to  HER2 
inactive  model  using  GOED.  The  ATP  binding  site  was  defined  as  the  active  site. 
Inspection  of  the  binding  modes  of  the  two  ligands  inside  the  ATP  binding  site  of  HER2 
models  showed  general  agreement  with  the  kinase  family  hinge  region  interactions  and 
correlates  with  established  pharmacophore  model.  To  analyze  the  binding  interactions  of 
molecules  with  HER2  active  and  inactive  states,  known  inhibitors  were  compared  with 
crystal  structures  of  HERl  active  and  inactive  states  (Eigures  7  and  8).  In  the  HER2 
active  state  complexed  with  erlotinib,  the  N(l)  of  quinazoline  ring  accepts  a  N-H---N 
bond  from  the  backbone  N-H  of  MetSOl  (Eigure  7).  A  strong  C-H---0  hydrogen  bond  is 
formed  between  the  activated  C(2)-H  group  and  the  backbone  carbonyl  oxygen  of 
Glu799.  The  distance  between  the  N3  atom  of  the  quinoline  and  OH  of  Thr798  is  3.54  A. 
Additionally,  the  6-  and  7-substituents  project  outside  the  active  site  (actually,  in  the  cleft 
region  between  the  C-terminal  and  N-terminal  lobes).  The  meta- substituted  4-anilino 
group  fits  well  in  the  hydrophobic  pocket  surrounded  by  amino  acid  residues  Phe864, 


Asp863,  Val797,  Leu796,  Met774,  Ile752,  Lys753,  Phe731.  In  the  complex  of  HER2 
inactive  state  and  Tykerb/GW2016  that  the  N1  atom  of  the  quinoline  is  involved  in  a  H- 

o 

bond  interaction  with  the  backbone  NH  of  MetSOl  at  a  distance  of  2.04  A.  The  distance 

o 

between  the  N3  atom  of  the  quinoline  and  OH  of  Thr798  is  3.39  A.  The  aniline  portion  of 
the  inhibitor  is  surrounded  by  residues  Thr798,  Lys753,  Thr862,  and  Asp841.  The 
pyridyl-methoxy  group  lies  in  a  pocket  surrounded  by  Val797,  Met774,  Arg784,  Ser783, 
Phe864,  Gly865  and  Leu866  (Figure  8).  It  was  observed  from  the  docking  studies  that  the 
overall  interactions  and  the  number  of  (C-H---0),  (0-H--N)  and  (N-H--N)  interactions 
of  ligands  that  are  forming  with  the  two  states  of  HER2  are  nearly  same.  A  second  type 
of  analysis  was  carried  out  to  verify  goodness  of  these  two  models  for  our  virtual 
screening  study  and  to  further  discriminate  the  performance  GOED,  Glide  and  eHITs 
programs.  A  diverse  small  set  of  molecules  (3,  7  -27)  was  docked  on  to  the  active  and 
inactive  states  of  HER2.  The  individual  scores  of  GOED,  Glide  and  eHITs  of  these 
twenty  two  compounds  onto  the  HER2  active  and  inactive  states  were  correlated  with 
experimental  activities  (pICjo)  (Table  4).  The  docking  results  and  the  correlation  of  their 
scores  with  experimental  activity  further  validated  our  models. 

Next,  we  wanted  to  know  which  score  correlates  better  with  experimental  activity 
of  HER2  antagonists.  In  order  to  answer  this  question,  the  docking  study  with  a  database 
of  474  known  HER2  inhibitors  on  HER2  active  and  inactive  states  were  carried  out.  The 
individual  scores  of  GOED,  Glide  and  eHITs  were  then  used  to  correlate  with 
experimental  pICjo  values  (Table  5).  Overall  we  found  that  docking  with  GOED,  Glide 
and  eHITs  performed  equally  well.  We  also  observed  that  docking  study  could  limit  false 


positives  but  showed  more  false  negatives.  To  overcome  the  false  positives  and  false 
negatives,  we  developed  a  combined  model  referred  as  a  consensus  model. 

Consensus  Model.  As  discussed  above,  no  single  scoring  function  and  parmacophore 
prediction  performed  accurately  enough  for  the  474  inhibitors  for  both  active  and  inactive 
state  of  HER2.  In  principle,  it  should  be  possible  to  determine  before  hand  which  score  is 
best  suited  in  a  given  case  and  then  to  use  this  score  to  decide  whether  a  docking  result  is 
likely  to  represent  a  hit  or  not.  However,  the  complexity  of  the  relationship  between  the 
performance  of  a  given  score  and  uniqueness  of  a  given  binding  site  and  the  docked 
compound  complicates  formulating  explicit  rules  for  such  a  relationship.  Therefore,  we 
chose  to  build  a  consensus  model  by  taking  predicted  pharmacophore  activity  and 
docking  scores  from  GOLD,  Glide  and  eHITs  obtained  from  active  and  inactive  states  of 
HER2  as  descriptors.^"^ 

The  final  equations,  Eq.  1  and  Eq.  2  were  generated  for  inactive  and  active  states  of 
HER2  using  multiple  linear  regressions 

Pred_activity  (pICjo)  =  0.233  +  0.7104  x  "Pharma_pIC5o"  -  0.063  x  "Glide" 

+  0.013  X  "Gold"-  0.159  x  "eHITS"  (1) 
n  =  158,  r^  =  0.846,  q^  =  0.684,  E  =  90.12,  PRESS  =  32.691  and  r^ed  =  0.69 

Pred_Activity  (pICjo)  =  0.206  -i-  0.698  x  "Pharma_pIC5o"  -  0.096  x  "Glide" 

-t  0.012  X  "Gold"-  0.127  x  "eHITS"  (2) 
n  =  158,  r^  =  0.763,  q^  =  0.674,  E  =  88.50,  PRESS  =  33.691  and  r^ed  =  0.66. 


Eq.  1  confirms  the  predicted  power  of  consensus  model  obtained  with  docking  scores  of 
HER2  inactive  state  along  with  pharmacophore  predicted  activity  as  proven  by 
significant  statistical  parameters.  Interestingly,  Eq.  2  obtained  with  docking  scores  of 
HER2  active  along  with  pharmacophore  predicted  activity  afforded  a  very  similar  and 
good  correlation  with  biological  activity  values  (Eigure  9).  The  scatter  plot  of  Eigure  11 
shows  remarkable  correlation  between  predicted  activity  values  from  consensus  models 
of  both  states  versus  experimental  biological  activity.  Eurther,  Eq.  1  and  Eq.  2  were  used 
to  calculate  a  number  of  statistical  parameters  such  as  hit  list  (Ht),  percent  active  (%Y), 
percent  ratio  of  actives  in  the  hit  list  (%A),  enrichment  factor  of  (E),  false  negatives,  false 
positives,  and  goodness  of  hit  score  to  cross  validate  the  predicted  power  of  the  models 
(Table  6). 

We  noticed  that  the  consensus  models  significantly  reduced  the  false  positives 
and  false  negatives  (Eigure  9  and  Table  6).  Equations  Eq.  1  and  Eq.  2  show  similar 
statistical  parameters  and  suggest  that  a  good  HER2  inhibitor  would  have  to  fit  with  the 
pharmacophore  and  have  to  successfully  interact  with  both  active  and  inactive  states  of 
HER2.  Therefore,  we  used  these  two  equations  to  screen  the  compounds  against  HER2. 

Compound  Selection  and  Cell  Growth  Inhibition  Assay.  The  compound  selection 
process  that  was  implemented  in  the  present  study  is  shown  as  a  flowchart  in  Figure  10. 
Einally,  57  compounds  were  selected  for  cytotoxicity  assay  in  SKBr3  cell  lines 
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overexpressing  HER2.  '  Tweleve  compounds  were  cytotoxic  with  IC50  values  ranging 
from  2  to  50  pM  and  are  presented  in  this  report  (Table  7). 


It  is  observed  from  the  Table  7  that  the  tricyclic  4-anilino  benzo[g] quinolines 
analogues  (1  and  2)  were  found  to  be  highly  potent  compounds  among  the  tested 
compounds.  The  tricyclic  4-anilino  benzo[g] quinoline  compound  1  and  2  inhibited  the 
proliferation  of  SKBr3  cell  lines,  which  highly  overexpresses  HER2,  with  IC50  value  of 
2.2  and  2.5pM,  respectively.  Compound  1  and  2  provide  activity  equivalent  to  that  of  the 
corresponding  the  tricyclic  3-cyanobenzo[g]quinoline  with  a  4-phenoxyanilino  group  at 
the  4-position,  which  were  earlier  reported  as  competitive  inhibitors  with  ATP  on  a  panel 
of  kinase  targets.  ’  It  is  also  observed  that  the  potency  of  the  compounds  decreases  with 
increasing  aliphatic  chain  between  benzo[g]quinolin-4-amine  and  phenyl  group  (3).  The 
replacement  of  phenyl  group  of  compound  3  with  a  1,4-methylpiperazine  led  to 
compound  4  with  decreased  inhibitory  activity.  As  envisaged  by  the  pharmacophore 
model,  most  of  the  2-methylquinoline  or  2-methylquinazoline  (5  -  6)  exhibited  enhanced 
inhibitory  activity  as  compared  to  compound  1  against  the  proliferation  of  SKBr3  cell 
lines. Partial  mapping  of  the  best  pharmacophore  model  onto  one  of  the  tested 
compounds,  a  quinoline  analogue  (7),  illustrates  that  two  carbon  linker  between  N  and 
phenyl  group  of  aniline  is  not  required  to  enhance  the  activity  of  the  compound. 
Compound  8,  containing  4-anilino  thieno[2,3-d]pyrimidin  core,  showed  with  IC50  values 
of  9  pM.  The  two  carbon  linker  containing  thieno[2,3-d]pyrimidin  analogue  (compound 
9)  inhibited  the  SKBr3  cell  growth  with  IC50  value  14.10  pM,  slightly  higher  than  that  of 
compound  8  (9.0  pM).  Surprisingly,  compound  10,  thieno[2,3-d]pyrimidine-4-thiol 
analogue,  showed  moderate  activity  with  IC50  value  of  16.5  pM.  The  [l,2,4]triazolo[l,5- 
a]pyrimidin-7-amine  scaffold  containing  compounds  (11  and  12)  displayed  moderate 


cytotoxic  activity  on  the  SKBr3  cell  lines  and  the  growth  of  the  cell  lines  was  inhibited 


with  IC50  values  of  9.0  and  4.8  |a,M,  respectively.  These  results  suggest  that  compound  1 
with  different  C-6  and  C-7  substituents  on  the  tricyclic  moiety  warrants  further 
modifications.  Further  structural  optimization  to  enhance  activity  is  in  progress  and  the 
results  will  be  described  in  the  due  course. 

Conclusions 

A  set  of  novel  cytotoxic  compounds  with  diverse  structural  scaffolds  were  identified 
using  two  consensus  models  Eq.  1  and  Eq.  2  that  were  generated  on  the  basis  of  known 
HER2  inhibitors.  Many  of  these  compounds  possess  amenable  chemical  and  structural 
features  are  useful  potential  leads  for  drug  design  strategies  targeting  HER2.  Eurther 
structural  optimization  to  enhance  HER2  inhibitory  activity  is  in  progress  and  will  be 
reported  in  due  course. 

Acknowledgment 

This  study  was  supported  in  part  by  DOD  Concept  Award  to  Nouri  Neamati.  Rambabu 
Gundla  would  like  to  thank  Sunita  Tajne,  Jinxia  Deng  and  Bikas  Debnath  for  their 
technical  support. 


References 

1.  Schlessinger,  J.  Cell  signaling  by  receptor  tyrosine  kinases.  Cell  2000,  103,  211- 
25. 

2.  Wells,  A.  EGE  receptor.  Int  J  Biochem  Cell  Biol  1999,  31,  637-43. 


3.  Guy,  P.  M.;  Platko,  J.  V.;  Cantley,  L.  C.;  Cerione,  R.  A.;  Carraway,  K.  L.,  3rd. 
Insect  cell-expressed  pl80erbB3  possesses  an  impaired  tyrosine  kinase  activity.  Proc 
Natl  Acad  Sci  U  S  A  1994,  91,  8132-6. 

4.  Jones,  J.  T.;  Akita,  R.  W.;  Sliwkowski,  M.  X.  Binding  specificities  and  affinities 
of  egf  domains  for  ErbB  receptors.  FEES  Lett  1999,  447,  227-31. 

5.  Cho,  H.  S.;  Mason,  K.;  Ramyar,  K.  X.;  Stanley,  A.  M.;  Gabelli,  S.  B.;  Denney,  D. 
W.,  Jr.;  Leahy,  D.  J.  Structure  of  the  extracellular  region  of  HER2  alone  and  in  complex 
with  the  Herceptin  Eab.  Nature  2003,  421,  756-60. 

6.  Garrett,  T.  P.;  McKern,  N.  M.;  Lou,  M.;  Elleman,  T.  C.;  Adams,  T.  E.;  Lovrecz, 
G.  O.;  Kofler,  M.;  Jorissen,  R.  N.;  Nice,  E.  C.;  Burgess,  A.  W.;  Ward,  C.  W.  The  crystal 
structure  of  a  truncated  ErbB2  ectodomain  reveals  an  active  conformation,  poised  to 
interact  with  other  ErbB  receptors.  Mol  Cell  2003,  11,  495-505. 

7.  Hynes,  N.  E.;  Lane,  H.  A.  ERBB  receptors  and  cancer:  the  complexity  of  targeted 
inhibitors.  Nat  Rev  Cancer  2005,  5,  341-54. 

8.  Moy,  B.;  Kirkpatrick,  P.;  Kar,  S.;  Goss,  P.  Lapatinib.  Nat  Rev  Drug  Discov  2007, 
6,  431-2. 

9.  Reid,  A.;  Vidal,  L.;  Shaw,  H.;  de  Bono,  J.  Dual  inhibition  of  ErbBl 
(EGER/HERl)  and  ErbB2  (HER2/neu).  Eur  J  Cancer  2007,  43,  481-9. 

10.  Spector,  N.;  Xia,  W.;  El-Hariry,  L;  Yarden,  Y.;  Bacus,  S.  HER2  therapy.  Small 
molecule  HER-2  tyrosine  kinase  inhibitors.  Breast  Cancer  Res  2007,  9,  205. 

11.  Moro,  S.;  Bacilieri,  M.;  Deflorian,  E.  Combining  ligand-based  and  structure-based 
drug  design  in  the  virtual  screening  arena.  Expert  Opin.  Drug  Discov.  2007,  2,  37-49. 

12.  Eeld,  R.;  Sridhar,  S.  S.;  Shepherd,  E.  A.;  Mackay,  J.  A.;  Evans,  W.  K.  Use  of  the 
epidermal  growth  factor  receptor  inhibitors  gefitinib  and  erlotinib  in  the  treatment  of  non¬ 
small  cell  lung  cancer:  a  systematic  review.  J  Thorac  Oncol  2006,  1,  367-76. 

13.  Kinase  Database.  Informatics,  GVK  Biosciences  Pvt.  Ltd,  S-1,  Phase- 1,  T.I.E. 
Balanagar,  Hyderabad  -  500  037,  India. 

14.  Gazit,  A.;  Osherov,  N.;  Posner,  L;  Yaish,  P.;  Poradosu,  E.;  Gilon,  C.;  Levitzki,  A. 
Tyrphostins.  2.  Heterocyclic  and  alpha- substituted  benzylidenemalononitrile  tyrphostins 
as  potent  inhibitors  of  EGE  receptor  and  ErbB2/neu  tyrosine  kinases.  J  Med  Chem  1991, 
34,  1896-907. 

15.  Singh,  J.;  Dobrusin,  E.  M.;  Ery,  D.  W.;  Haske,  T.;  Whitty,  A.;  McNamara,  D.  J. 
Structure-based  design  of  a  potent,  selective,  and  irreversible  inhibitor  of  the  catalytic 
domain  of  the  erbB  receptor  subfamily  of  protein  tyrosine  kinases.  J  Med  Chem  1997,  40, 
1130-5. 

16.  Sun,  L.;  Tran,  N.;  Tang,  E.;  App,  H.;  Hirth,  P.;  McMahon,  G.;  Tang,  C.  Synthesis 
and  biological  evaluations  of  3-substituted  indolin-2-ones:  a  novel  class  of  tyrosine 
kinase  inhibitors  that  exhibit  selectivity  toward  particular  receptor  tyrosine  kinases.  J 
Med  Chem  1998,  41,  2588-603. 

17.  Smaill,  J.  B.;  Rewcastle,  G.  W.;  Loo,  J.  A.;  Greis,  K.  D.;  Chan,  O.  H.;  Reyner,  E. 
L.;  Lipka,  E.;  Showalter,  H.  D.;  Vincent,  P.  W.;  Elliott,  W.  L.;  Denny,  W.  A.  Tyrosine 
kinase  inhibitors.  17.  Irreversible  inhibitors  of  the  epidermal  growth  factor  receptor:  4- 
(phenylamino)quinazoline-  and  4-(phenylamino)pyrido[3,2-d]pyrimidine-6-acrylamides 
bearing  additional  solubilizing  functions.  J  Med  Chem  2000,  43,  1380-97. 


18.  Wells,  G.;  Seaton,  A.;  Stevens,  M.  F.  Structural  studies  on  bioactive  compounds. 
32.  Oxidation  of  tyrphostin  protein  tyrosine  kinase  inhibitors  with  hypervalent  iodine 
reagents.  J  Med  Chem  2000,  43,  1550-62. 

19.  Tsou,  H.  R.;  Mamuya,  N.;  Johnson,  B.  D.;  Reich,  M.  F.;  Gruber,  B.  C.;  Ye,  F.; 
Nilakantan,  R.;  Shen,  R.;  Discafani,  C.;  DeBlanc,  R.;  Davis,  R.;  Koehn,  F.  E.; 
Greenberger,  L.  M.;  Wang,  Y.  F.;  Wissner,  A.  6-Substituted-4-(3- 
bromophenylamino)quinazolines  as  putative  irreversible  inhibitors  of  the  epidermal 
growth  factor  receptor  (EGFR)  and  human  epidermal  growth  factor  receptor  (HER- 2) 
tyrosine  kinases  with  enhanced  antitumor  activity.  J  Med  Chem  2001,  44,  2719-34. 

20.  Smaill,  J.  B.;  Showalter,  H.  D.;  Zhou,  H.;  Bridges,  A.  J.;  McNamara,  D.  J.;  Ery, 
D.  W.;  Nelson,  J.  M.;  Sherwood,  V.;  Vincent,  P.  W.;  Roberts,  B.  J.;  Elliott,  W.  E.; 
Denny,  W.  A.  Tyrosine  kinase  inhibitors.  18.  6-Substituted  4-anilinoquinazo lines  and  4- 
anilinopyrido[3,4-d]pyrimidines  as  soluble,  irreversible  inhibitors  of  the  epidermal 
growth  factor  receptor.  J  Med  Chem  2001,  44,  429-40. 

21.  Wissner,  A.;  Overbeek,  E.;  Reich,  M.  E.;  Eloyd,  M.  B.;  Johnson,  B.  D.;  Mamuya, 
N.;  Rosfjord,  E.  C.;  Discafani,  C.;  Davis,  R.;  Shi,  X.;  Rabindran,  S.  K.;  Gruber,  B.  C.; 
Ye,  E.;  Hallett,  W.  A.;  Nilakantan,  R.;  Shen,  R.;  Wang,  Y.  E.;  Greenberger,  E.  M.;  Tsou, 
H.  R.  Synthesis  and  structure- activity  relationships  of  6,7-disubstituted  4- 
anilinoquinoline-3-carbonitriles.  The  design  of  an  orally  active,  irreversible  inhibitor  of 
the  tyrosine  kinase  activity  of  the  epidermal  growth  factor  receptor  (EGER)  and  the 
human  epidermal  growth  factor  receptor-2  (HER- 2).  J  Med  Chem  2003,  46,  49-63. 

22.  Tsou,  H.  R.;  Overbeek- Klumpers,  E.  G.;  Hallett,  W.  A.;  Reich,  M.  E.;  Eloyd,  M. 
B.;  Johnson,  B.  D.;  Michalak,  R.  S.;  Nilakantan,  R.;  Discafani,  C.;  Golas,  J.;  Rabindran, 
S.  K.;  Shen,  R.;  Shi,  X.;  Wang,  Y.  E.;  Upeslacis,  J.;  Wissner,  A.  Optimization  of  6,7- 
disubstituted-4-(arylamino)quinoline-3-carbonitriles  as  orally  active,  irreversible 
inhibitors  of  human  epidermal  growth  factor  receptor-2  kinase  activity.  J  Med  Chem 
2005,  48,  1107-31. 

23.  Elauger,  E.;  He,  H.;  Kim,  J.;  Aguirre,  J.;  Rosen,  N.;  Peters,  U.;  Davies,  P.; 
Chiosis,  G.  Evaluation  of  8-arylsulfanyl,  8-arylsulfoxyl,  and  8-arylsulfonyl  adenine 
derivatives  as  inhibitors  of  the  heat  shock  protein  90.  J  Med  Chem  2005,  48,  2892-905. 

24.  Borzilleri,  R.  M.;  Zheng,  X.;  Qian,  E.;  Ellis,  C.;  Cai,  Z.  W.;  Wautlet,  B.  S.; 
Mortillo,  S.;  Jeyaseelan,  R.,  Sr.;  Kukral,  D.  W.;  Eura,  A.;  Kamath,  A.;  Vyas,  V.; 
Tokarski,  J.  S.;  Barrish,  J.  C.;  Hunt,  J.  T.;  Eombardo,  E.  J.;  Eargnoli,  J.;  Bhide,  R.  S. 
Design,  synthesis,  and  evaluation  of  orally  active  4-(2,4-difluoro-5- 
(methoxycarbamoyl)phenylamino)pyrrolo[2,l-f][l,2,4]triaz  ines  as  dual  vascular 
endothelial  growth  factor  receptor-2  and  fibroblast  growth  factor  receptor- 1  inhibitors.  J 
Med  Chem  2005,  48,  3991-4008. 

25.  Gaul,  M.  D.;  Guo,  Y.;  Affleck,  K.;  Cockerill,  G.  S.;  Gilmer,  T.  M.;  Griffin,  R.  J.; 
Guntrip,  S.;  Keith,  B.  R.;  Knight,  W.  B.;  Mullin,  R.  J.;  Murray,  D.  M.;  Rusnak,  D.  W.; 
Smith,  K.;  Tadepalli,  S.;  Wood,  E.  R.;  Eackey,  K.  Discovery  and  biological  evaluation  of 
potent  dual  ErbB-2/EGER  tyrosine  kinase  inhibitors:  6-thiazolylquinazolines.  Bioorg 
Med  Chem  Lett  2003,  13,  637-40. 

26.  Revesz,  E.;  Blum,  E.;  Di  Padova,  E.  E.;  Buhl,  T.;  Eeifel,  R.;  Gram,  H.;  Hiestand, 
P.;  Manning,  U.;  Rucklin,  G.  Novel  p38  inhibitors  with  potent  oral  efficacy  in  several 
models  of  rheumatoid  arthritis.  Bioorg  Med  Chem  Lett  2004,  14,  3595-9. 


27.  Bridges,  A.  J.  Chemical  inhibitors  of  protein  kinases.  Chem  Rev  2001,  101,  2541- 
72. 

28.  Garcia-Echeverria,  C.;  Fabbro,  D.  Therapeutically  targeted  anticancer  agents: 
inhibitors  of  receptor  tyrosine  kinases.  Mini  Rev  Med  Chem  2004,  4,  273-83. 

29.  Stuart,  C.  G.;  Clive,  C.  M.;  Karl,  M.  S.;  Sadie,  V.;  John,  P.  M.;  Thomas,  H.  A.; 
Paul,  B.;  Witold,  F.  K.  Preparation  of  heterocyclyl- substituted  quinazo lines  as  protein 
tyrosine  kinase  inhibitors.  WO  9703069  Al,  1997. 

30.  Stuart,  C.  G.;  Clive,  C.  M.;  Barry,  G.  S.;  Jane,  S.  K.  Preparation  of  bicyclic 
heteroaromatic  compounds  as  protein  tyrosine  kinase  inhibitors.  WO  9802437  Al,  1998. 

31.  Ronghui,  F.;  Peter,  J.  C.;  Steven,  W.;  Shenlin,  H.;  Stuart,  E.;  Robert,  G.;  Steve, 
M.  Preparation  of  l,2,4-triazole-3, 5-diamine  derivatives  as  kinase  inhibitors.  WO 
2002057240  Al,  2002. 

32.  Peter,  H.  K.;  Pruess,  S.  D.;  Elaina,  M.;  Kay,  S.  F.;  Gyorgi,  K.;  Istvan,  S.;  Tamas, 
B.;  Janis,  H.;  Faszlo,  O.;  Alex,  F.;  Aviv,  G.;  Axel,  U.;  Reiner,  F.;  Fairooz,  F.  K.;  Dennis, 
S.;  Cho,  T.  P.  Treatment  of  platelet  derived  growth  factor-related  disorders  such  as 
cancers.  US  5958959  A,  1999. 

33.  Chen,  H.;  Gazit,  A.;  Fevitzki,  A.;  Hirth,  K.  P.;  Mann,  E.;  Shawver,  K.  F.;  Tsai,  J.; 
Tang,  P.  C.  Methods  and  compositions  for  inhibiting  cell  proliferative  disorders.  US 
20020068687  Al,  2002. 

34.  Yu,  M.;  Etsuya,  M.  Preparation  of  aralkylazoles  as  tyrosine  kinase  inhibitors 
useful  as  antitumor  agents.  WO  9803505  A2,  1998. 

35.  Akihiro,  T.;  Takenori,  H.;  Etsuya,  M.  Preparation  of  2-styryl-4- 
(phenoxymethyl)oxazole  derivatives  as  tyrosine  kinase  inhibitors.  WO  2001077107  Al, 

2001. 

36.  Hui,  C.;  Aviv,  G.;  Peter,  K.  H.;  Elaina,  M.;  Faura,  K.  S.;  Jianming,  T.;  Cho,  P.  T. 
Methods  and  compositions  using  receptor  tyrosine  kinase  inhibitors  for  inhibiting  cell 
proliferative  disorders,  and  inhibitor  preparation.  US  5789427  A,  1998. 

37.  Neal,  R.;  Scott,  D.  K.;  Samuel,  J.  D.;  Furzhong,  F.  Z.;  Faura,  S.;  Ouathek,  O. 
Methods  and  compositions  using  bifunctional  hsp-binding  derivatives  for  degradation 
and/or  inhibition  of  HER-family  tyrosine  kinases  and  treatment  of  cancer.  WO 
2000061578  Al,  2000. 

38.  Mien-Chie,  H.;  Fisha,  Z.  Emodin  and  related  compounds  for  sensitization  of 
HER2/neu  over-expressing  cancer  cells  to  chemotherapeutic  drugs.  WO  9727848  Al, 

1997. 

39.  Stuart,  G.  C.;  Clive,  M.  C.;  Barry,  S.  G.;  Jane,  K.  S.  Preparation  of 
azolylquinazo lines  and  related  compounds  as  protein  tyrosine  kinase  inhibitors.  WO 
9802434  Al,  1998. 

40.  Thomas,  A.  H.;  Sadie,  V.;  Paul,  B.;  Karl  Witold,  F.;  Stephen  Carl,  M.;  P.,  M.  J. 
Preparation  of  quinoline  and  quinazo  line  protein  tyrosine  kinase  inhibitors.  WO  9609294 

Al,  1996. 

41.  Hugh,  R.  B.;  Grant,  J.  K.;  Stewart,  J.  S.;  Christophe,  B.  B.  Preparation  of  4- 
anilino-5-carbamoylmethoxyquinazolines  as  selective  erbB2  receptor  tyrosine  kinase 
inhibitors  for  use  against  tumors.  WO  2005118572  Al,  2005. 

42.  Cerius2.  4.11;  Accelrys  Inc:  San  Diego,  CA,  USA,  2005. 

43.  Catalyst.  4.11;  Accelrys  Inc:  San  Diego,  CA,  USA,  2005. 


44.  Smellie,  A.  T.,  S.  L.;  Towbin,  P.  Polings-promoting  conformational  variation.  J. 
Comput.  Chem  1995,  16,  171-187. 

45.  Stamos,  J.;  Sliwkowski,  M.  X.;  Eigenbrot,  C.  Structure  of  the  epidermal  growth 
factor  receptor  kinase  domain  alone  and  in  complex  with  a  4-anilinoquinazoline  inhibitor. 
J  Biol  Chem  2002,  277,  46265-72. 

46.  Wood,  E.  R.;  Truesdale,  A.  T.;  McDonald,  O.  B.;  Yuan,  D.;  Hassell,  A.; 
Dickerson,  S.  H.;  Ellis,  B.;  Pennisi,  C.;  Horne,  E.;  Eackey,  K.;  Alligood,  K.  J.;  Rusnak, 
D.  W.;  Gilmer,  T.  M.;  Shewchuk,  E.  A  unique  structure  for  epidermal  growth  factor 
receptor  bound  to  GW572016  (Eapatinib):  relationships  among  protein  conformation, 
inhibitor  off-rate,  and  receptor  activity  in  tumor  cells.  Cancer  Res  2004,  64,  6652-9. 

47.  Eiao,  J.  J.  Molecular  recognition  of  protein  kinase  binding  pockets  for  design  of 
potent  and  selective  kinase  inhibitors.  J  Med  Chem  2007,  50,  409-24. 

48.  GOED.  1.2;  CCDC:  Cambridge,  UK,  2002. 

49.  Eriesner,  R.  A.;  Murphy,  R.  B.;  Repasky,  M.  P.;  Erye,  E.  E.;  Greenwood,  J.  R.; 
Halgren,  T.  A.;  Sanschagrin,  P.  C.;  Mainz,  D.  T.  Extra  precision  glide:  docking  and 
scoring  incorporating  a  model  of  hydrophobic  enclosure  for  protein- ligand  complexes.  J 
Med  Chem  2006,  49,  6177-96. 

50.  Zsoldos,  Z.;  Reid,  D.;  Simon,  A.;  Sadjad,  S.  B.;  Johnson,  A.  P.  eHiTS:  A  new 
fast,  exhaustive  flexible  ligand  docking  system.  J  Mol  Graph  Model  2007,  26,  198-212. 

51.  Eantin,  V.  R.;  Berardi,  M.  J.;  Babbe,  H.;  Michelman,  M.  V.;  Manning,  C.  M.; 
Eeder,  P.  A  bifunctional  targeted  peptide  that  blocks  HER-2  tyrosine  kinase  and  disables 
mitochondrial  function  in  HER-2-positive  carcinoma  cells.  Cancer  Res  2005,  65,  6891- 
900. 

52.  Gril,  B.;  Vidal,  M.;  Assayag,  E.;  Poupon,  M.  E.;  Eiu,  W.  Q.;  Garbay,  C.  Grb2- 
SH3  ligand  inhibits  the  growth  of  HER2-I-  cancer  cells  and  has  antitumor  effects  in  human 
cancer  xenografts  alone  and  in  combination  with  docetaxel.  Int  J  Cancer  2007,  121,  407- 
15. 

53.  Jin,  Y.;  Ei,  H.  Y.;  Ein,  E.  P.;  Tan,  J.;  Ding,  J.;  Euo,  X.;  Eong,  Y.  Q.  Synthesis  and 
antitumor  evaluation  of  novel  5-substituted-4-hydroxy-8-nitroquinazolines  as  EGER 
signaling-targeted  inhibitors.  Bioorg  Med  Chem  2005,  13,  5613-22. 

54.  Aparna,  V.;  Rambabu,  G.;  Panigrahi,  S.  K.;  Sarma,  J.  A.;  Desiraju,  G.  R.  Virtual 
screening  of  4-anilinoquinazoline  analogues  as  EGER  kinase  inhibitors:  importance  of 
hydrogen  bonds  in  the  evaluation  of  poses  and  scoring  functions.  J  Chem  Inf  Model  2005, 
45,  725-38. 

55.  Zhang,  N.;  Wu,  B.;  Wissner,  A.;  Powell,  D.  W.;  Rabindran,  S.  K.;  Kohler,  C.; 
Boschelli,  E.  4-Anilino-3-cyanobenzo[g]quinolines  as  kinase  inhibitors.  Bioorg  Med 
Chem  Lett  2002,  12,  423-5. 


Figure  legends 

Figure  1:  Clinically  used  HER  family  targeted  kinase  inhibitors. 

Figure  2:  Training  set  compounds  (4  -  9)  used  to  generate  qualitative  pharmacophore 
models. 

Figure  3:  Training  set  compounds  (10  -  27)  used  to  generate  quantitative  pharmacophore 
models. 

Figure4:  ClustalW  alignment  of  HERl  (EGER)  and  HER2  with  highlighted  residues  of 
Asp-746,  Cys-751  and  DEG  motif. 

Figure  5:  A)  The  three  dimensional  arrangement  of  pharmacophoric  features  in  the 

quantitative  pharmacophore  model  (Hypo5).  Green:  hydrogen  bond  acceptor 
(HBA),  blue:  hydrophobic  aromatics  (HRAl-  HRA3),  magenta:  hydrogen 
bond  donor  (HBD).  B)  Mapping  of  Hypo5  onto  Tykerb  (3).  C)  Mapping  of 
Hypo5  onto  compound  14  from  the  moderately  active  category,  D)  Mapping  of 
Hypo5  onto  compound  20  from  the  minimally  active  category. 

Figure  6:  Scatter  plot  of  experimental  vs  predicted  activities  using  pharmacophore  model 
Hypo5. 

Figure  7:  A)  HERl  active  state  conformation  with  Erlotinib,  B)  HER2  active  state 
conformation  with  Erlotinib. 

Figure  8:  A)  HERl  inactive  state  conformation  with  Tykerb/GW2016,  B)  HER2  inactive 
state  conformation  with  Tykerb/GW2016. 

Figure  9:  A)  Consensus  model  of  HER2  inactive  conformation  with  known  HER2 
inhibitors,  B)  Consensus  model  of  HER2  active  conformation  with  known 


HER2  inhibitors. 


Figure  10:  In  silico  screening  protocol  implemented  in  the  discovery  of  HER2  inhibitors. 
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Tablel.  Experimental  and  Predicted  IC50  Values  of  22  Training  Set  compounds  with 
Hypo5 


Compound 

Experimental 
IC50  (nM) 

Predicted 

IC50  (nM) 

Pit  value 

7 

1 

1.5 

9.49 

8 

1 

1.4 

9.51 

10 

2 

3 

9.20 

11 

4 

7.5 

8.79 

3 

9 

5.7 

8.91 

12 

10 

8.7 

8.73 

13 

13 

47 

7.99 

14 

13 

14 

8.51 

9 

19 

6.7 

8.84 

15 

120 

270 

7.24 

16 

200 

180 

7.42 

17 

733 

190 

7.28 

18 

1000 

380 

7.09 

19 

1700 

3990 

6.69 

20 

5300 

1800 

6.40 

21 

9300 

3900 

6.08 

22 

16000 

38000 

5.09 

23 

21000 

89000 

4.72 

24 

35000 

47000 

5.00 

25 

65000 

53000 

4.95 

26 

120000 

71000 

4.82 

27 

310000 

40000 

5.06 

Table  2.  Statistical  Parameters  of  Ten  Pharmacophore  Models  Generated  by  the  Hypogen 
Algorithm 


Hypo 

No. 

Total  cost 

Cost 

difference^* 

Error 

Cost 

RMS 

deviation 

Training  set 

(r) 

Eeatures'’ 

1 

78.639 

80.301 

62.598 

0.798 

0.973 

ADRR 

2 

79.722 

79.218 

63.796 

0.882 

0.967 

ADRR 

3 

81.577 

77.362 

65.254 

0.974 

0.960 

ADRR 

4 

81.636 

77.303 

66.339 

1.037 

0.953 

ADRRR 

5 

81.791 

77.148 

66.456 

1.044 

0.953 

ADRRR 

6 

81.874 

77.065 

64.108 

0.902 

0.966 

ADRR 

7 

82.846 

76.093 

66.451 

1.044 

0.953 

ADRR 

8 

83.890 

75.049 

67.846 

1.120 

0.946 

ADRRR 

9 

85.148 

73.791 

69.858 

1.221 

0.935 

ADRRR 

10 

86.974 

71.965 

68.504 

1.154 

0.944 

ADRR 

‘“(Null  cost-Total  cost),  Null  cost  =  158.94,  Fixed  cost  =  72.47, 

For  Hypo-5  Weight  =  1.169,  Configuration  =  14.164. 

'’A-  Hydrogen  Bond  Acceptor,  D-  Hydrogen  Donor  and  R  -  Hydrophobic  aromatic. 


Table  3.  Hypo5  Pharmacophore  Model  Validation  with  a  Database  of  Known  HER2 
Inhibitors 


No 

Parameter 

Value 

1 

Total  compounds  in  database  (D) 

474 

2 

Total  number  of  actives  in  database  (A) 

195 

3 

Total  hits  (Ht) 

204 

4 

Active  hits  (Ha) 

140 

5 

%  Yield  of  actives 

68.63 

6 

%  Ratio  of  actives  in  the  hit  list 

71.79 

7 

Enrichment  factor  or  enhancement  (E) 

1.67 

8 

Ealse  negatives 

55 

9 

Ealse  positives 

64 

10 

GH  score  (Goodness  of  hit  list)** 

0.536 

‘“(Ha/4HtA)(3A-i-Ht))x(l-((Ht-Ha)/(D-A))];  GH  Score  of  0.6  -  0.7  indicates  a  very  good 
model. 


Table  4.  Correlation  Coefficient  Values  of  GOLD  score,  Glide  score,  cHITs  score 

(HER2  Active  and  Inactive  States)  with  Experimental  pICso  of  Known  22  HER2 
Inhibitors  (table  1) 


GOED 

Glide 

eHITS 

Inactive  state^* 

0.844 

-0.814 

-0.695 

Active  state^* 

0.805 

-0.704 

-0.790 

‘“HER2  conformation. 

Table  5.  Correlation  Coefficient  Values  of  GOED  Score,  Glide  score,  cHITs  score 
(HER2  Active  and  Inactive  States)  with  Experimental  pICso  of  Known  474 
HER2  Inhibitors 


GOED 

Glide 

CHITS 

Inactive  state^* 

0.551 

-0.596 

-0.549 

Active  state^* 

0.545 

-0.589 

-0.555 

‘“HER2  conformation. 

Table  6.  A)  Consensus  Model  of  HER2  Inactive  State  Conformation  and  Validation  with 
Known  HER2  Inhibitors  Database,  B)  Consensus  Model  of  HER2  Active 
State  Conformation  and  Validation  with  Known  HER2  Inhibitors  Database 


A) 

B) 

S.  No 

Parameter 

Value 

S.  No 

Parameter 

Value 

1 

Total  compounds  in  database  (D) 

474 

1 

Total  compounds  in  database  (D) 

474 

2 

Total  number  of  actives  in 
database  (A) 

195 

2 

Total  number  of  actives  in 
database  (A) 

195 

3 

Total  hits  (HO 

193 

3 

Total  hits  (HO 

197 

4 

Active  hits  (Ha) 

156 

4 

Active  hits  (Ha) 

154 

5 

%  Yield  of  actives 

81 

5 

%  Yield  of  actives 

78 

6 

%  Ratio  of  actives  in  the  hit  list 

80 

6 

%  Ratio  of  actives  in  the  hit  list 

79 

7 

Enrichment  or  Enhancement 
factor  (E) 

1.96 

7 

Enrichment  or  Enhancement 
factor  (E) 

1.90 

8 

Ealse  negatives 

39 

8 

Ealse  negatives 

41 

9 

Ealse  positives 

37 

9 

Ealse  positives 

43 

10 

GH  score  (Goodness  of  hit  list) 

0.699 

10 

GH  score  (Goodness  of  hit  list) 

0.663 

[(Ha/4HtA)(3A+Ht))x(l-((Ht-Ha)/(D-A))];  GH  Score  of  0.6  -  0.7  indicates  a  very  good 
model. 


Table  7.  Cytotoxicity  of  Novel  Compounds  Against  a  SKBr3  Breast  Cancer  Cells 
Overexpressing  HER2). 

Predicted  Experimental 

IC50  (pM)  IC50  (pM) 

S.No  Structure  _ 

Eql  Eq2  SKBr3  cell  lines 


0.71 


0.56 


16.9 


55.3 


0.47 


2.2 +  0.2 


2.5 +  0.2 


8.5 +  0.3 


15.02+  0.4 


5.1  +0.1 


6 


0.34 


0.53 


5  +  0.2 


7 


37.1 


41.17 


15.05  +  0.5 


0.51 


25.3 


24.12 


24.1 


49.6 


0.75 


32.7 


22.17 


26.3 


62.5 


9  +  0.5 


14.10  +  0.5 


16.5  +  0.3 


9  +  0.5 


14.08  +  0.4 


